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SUMMARY 


Measurements  have  been  made  of  the  instantaneous  and  time- 
averaged  values  of  the  induced  velocities  in  the  wake  of  a 
full-scale  helicopter  in  actual  steady  flight  conditions. 
These  measurements  were  made  with  three-dimensional  pitot 
static  and  two-dimensional  hot  film  probes.  The  test  setup 
provided  for  measurements  at  several  radial ,  axial,  and 
azimuthal  positions  relative  to  the  rotor  centerline. 

Simple,  analytical  methods  are  formulated  to  predict  the 
instantaneous  and  time -aver aged,  induced  velocities  in  the 
rotor  wake. 

Calculated  values  of  the  induced  velocities  are  compared  with 
the  experimental  data,  and  agreement  is  obtained  for  the 
hovering  flight  condition. 
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fuselage  equivalent  radius,  ft. 
a/R 

rotor  diameter,  ft. 

vertical  distance  from  rotor  plane  to  ground,  ft. 
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(see  Figure  1) 
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rotor  thrust,  lb. 

resultant  axial  velocity  through  the  rotor  disk, 
lb. /sec. 

instantaneous  induced  velocity  inside  the  rotor 
wake,  ft. /sec. 


u  time-averaged,  induced  velocity  inside  the 

rotor  wake,  ft. /sec. 

u0  momentum  value  of  rotor  induced  velocity  at 

the  rotor  disk,  ft. /sec. 

*  axial  distance  from  the  rotor  disk,  ft. 
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*  vertical  distance  from  the  plane  containing  the 
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(see  Figure  8) 

*'  x/R 

a  angular  position  of  any  point  P  with  respect 

to  the  positive  Y-axis  (see  Figure  1) 

a,  angle  defining  point  P  with  vertex  on  X^-axis 

parallel  to  the  Y^-Z^  plane 

a  angle  defining  point  P  with  vertex  on  X^-axis 

measured  from  the  positive  Y-axis  parallel  to 
the  X-Y  plane 

r  circulation  of  the  vortex  system  representing 

one  blade,  ft.  /sec. 

Tf  blade  tip  circulation  of  a  triangularly  loaded 

blade,  ft. 2/sec. 

9  angular  position  of  blade  with  respect  to  the 

Y-axis  (see  Figure  1) 

9  angular  position  of  any  point  P  with  respect 

to  the  blade, 9'=9-a  (see  Figure  1) 

p  air  mass  density,  slug/ft. 3 

t  time,  sec. 
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$  yawed  angle  (see  Figure  6) 

X  wake  angle  (see  Figure  6) 

^  azimuth  test  position  measured  from  the  downwind 

position  in  the  direction  of  cu 

stream  function  of  the  three-dimensional  doublet 
rep^  senting  the  fuselage 

*^n  stream  function  of  the  induced  flow  inside  the 

wake  of  a  free  lifting  rotor 

( v  rotor  angular  speed,  rad. /sec. 
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I.  INTRODUCTION 


The  induced  flow  is  one  of  the  important  parameters  which 
determine  the  aerodynamic  loads  of  the  fuselage  and  wing 
surfaces  immersed  in  the  wake  of  helicopter  rotors.  Large 
effects  on  the  stability  and  control  of  helicopters  can  be 
experienced  as  a  result  of  the  induced  flow  acting  on  tail 
surfaces.  Very  pronounced  effects  of  the  induced  flow  field 
are  also  evident  in  the  performance,  stability,  and  vibration 
characteristics  of  tandem  rotor  helicopters,  where  the  rear 
rotor  blades  operating  in  the  wake  of  the  front  rotor  down- 
wash  are  subjected  to  local  alternating  air  loads.  The 
actual  flow  field  in  the  wake  of  helicopter  rotors  varies  both 
with  time  and  location  relative  to  the  rotor.  Previous 
analyses  of  these  alternating  air  velocities  are  relatively 
complex  and  require  the  use  of  considerable  time  of  high-speed 
digital  computers.  Although  digital  computer  methods  repre¬ 
sent  an  extremely  useful  tool  for  some  applications, 
especially  since  they  permit  the  elimination  of  several 
restrictive  assumptions,  it  is  recognized  that  the  use  of 
these  computers  is  an  expensive  proposition.  It  is  desirable, 
therefore,  to  continue  the  search  for  analytical  methods  which 
promise  to  provide  a  simpler  and  preferably  a  closed-form  type 
of  solution  to  the  rotor  wake  velocity  problem. 

The  effort  under  the  program  reported  herein  was  directed, 
therefore,  toward  the  development  such  simpler  analytical 
methods.  In  addition,  full-scale  measurements  were  performed 
of  the  time-averaged  and  time-dependent  velocities  at  several 
locations  in  the  rotor  wake  of  a  Hughes  Model-269  helicopter. 

Comparisons  were  made  of  a  number  of  test  points  and  the 
corresponding  theoretical  values.  The  results  of  these 
comparisons  as  well  as  the  details  of  the  theoretical  analyses 
and  the  test  program  are  described  in  the  following  sections. 
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II.  THEORETICAL  ANALYSIS 


The  aerodynamic  forces  acting  on  the  fuselage  and  lifting 
surfaces  of  helicopters  are  greatly  influenced  by  the  rotor 
induced  flow.  This  flow,  in  general,  is  time  dependent  with 
frequencies  which  are  harmonics  of  the  rotor  angular  speed. 
Available  analytical  methods,  such  as  those  of  References  1 
through  3,  are  relatively  complex  and  require  the  use  of 
digital  computers  to  obtain  specific  solutions. 

In  Reference  4  an  analysis  is  presented  which  results  in 
simple,  closed-form  solutions  to  the  normal  average  induced 
velocity  in  the  wake  of  a  rotor  for  both  hover  and  forward 
speeds . 

In  the  program  reported  herein,  the  analysis  of  Reference  4 
has  been  extended  to  obtain  the  instantaneous,  in  addition 
to  the  time-averaged,  induced  velocity  in  the  rotor  wake. 
Also,  an  attempt  has  been  made  to  account  for  the  effect  of 
the  fuselage  on  the  induced  flow. 

The  analysis  is  performed  in  several  parts  as  follows: 

A.  Axial  Flow  D.  Ground  Effect 

B.  Yawed  Flow  E.  Wake  Contraction 

C.  Fuselage  Effects 

A.  AXIAL  FLOW 

The  average  velocity,  U  ,  through  a  rotor  or  a  propeller 
operating  in  an  axial  flow  condition  can  be  expressed  as 


U  =  V  +  u0 


(1) 


where  V  is  the  free-stream  velocity,  and  uq  is  the  momentum 
value  of  rotor  induced  velocity  at  the  rotor  disk,  given  by 


T _ 

2ptr  R2  (V  +u0  ) 


(2) 
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For  any  given  axial  flight  condition,  U  and  u0  can  therefore 
be  determined  from  the  above  two  equations. 

In  the  present  analysis,  the  normal  induced  velocity  in  the 
rotor  wake  is  derived  for  two  idealized  loadings,  i.e., 
uniform  and  triangular.  The  analysis  can,  however,  be 
extended  to  the  general  case  of  arbitrary  loading. 

The  velocity  induced  at  point  P  inside  the  rotor  wake  is 
derived  with  respect  to  a  cylindrical  coordinate  system 
fixed  to  the  n^  blade.  Definitions  of  this  coordinate 
system  and  the  induced  velocity  components  are  illustrated 
in  Figure  1.  As  noted  in  Figure  2,  the  effect  of  rotor  wake 
contraction  is  neglected. 

1 .  Uniform  Disk  Loading 

A  uniformly  loaded  rotor  is  also  one  of  constant  circulation 
along  the  rotor  blade  span.  In  this  case,  vorticity  is  shed 
only  at  the  blade  tips  and  the  rotor  axis.  The  induced 
velocities  are  derived  based  on  the  wake  structure  shown  in 
Figure  2.  According  to  this  wake  structure,  each  blade  is 
approximated  by  a  vortex  system  consisting  of  a  bound  vortex 
fixed  to  the  blade,  a  vortex  helix  which  trails  from  the 
blade,  and  a  central  vortex  along  the  rotor  axis.  This 
vortex  line  representation  is  similar  to  that  used  in  the 
theory  of  finite  wings  having  a  large  aspect  ratio 
(Reference  5). 

a.  Time-Averaged,  Induced  Velocity 

(i)  Normal  Velocity  Component,  un 

The  time-averaged,  induced  velocity  is  derived 
in  Reference  4  for  an  incompressible,  irrota- 
tional,  and  nonviscous  flow.  Using  Reference  4, 
we  find  that  the  perturbation  pressure  can  be 
expressed  as 


P,  =  -?u(uNh+UNb)"^u;r{uah+  uac*  uab) 


(3) 
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For  the  time-averaged  pressure  condition,  the 
effect  of  the  bound  vortex  can  be  neglected; 
hence,  equation  (3)  reduces  to 


p,=-/°UuNh“/°wr(uah  +  uac) 


When  the  notation  of  Figure  2  is  used,  the 
induced  velocity  at  any  point  P(x,y,z)  can  be 
computed  by  use  of  the  Biot-Savart  law 


ds  x  5 
la  I3 


(5) 


where  ds 
element , 
from  the 
velocity 


is  the  vector  length  of  a  vortex 
a  is  the  vector  distance  to  point  P 
vortex  element,  and  dq  is  the  induced 
vector. 


The  normal  velocity  component  is  induced  only  by 
the  helical  vortex.  As  shown  in  Reference  4, 
this  velocity  component  can  be  written  as 


r  ®[l-r'cos  (0-^)]de’ 

UNh=  4JR  J  }  ~2  I  7 

0  |(X  -  O2  + r'  + 1  -  2  r'  cos  (6-  )  j 


(6) 


3/2 


where 


Z'  = 

ut 

R 

e-a 

1  - 

u 

J  - 

nD 
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Similarly,  the  tangential  velocity  component 
induced  by  the  helical  vortex  is 


r 

Uah"47rR 


/ 


0 


[r'-tx'-Oysin  (0'-jO_cos 

It'Z'f 2r,cos(0,--je')  V2 


(7) 


Since  neither  equation  (6)  nor  equation  (7)  is 
integrable,  the  method  proceeds  by  first 
substituting  equations  (6)  and  (7)  into 
equation  (4).  There  follows  that 

p?u  ®  [l-2r‘  cos  (e'-^O+^-rV-O-J  sin (6  -  Jf'ldf 
4\  [(x'-Oz+r'%l-2r'cos  <e'-y{')]3/2 

-pur'  uac 

It  can  be  shown  by  direct  differentiation  that 
the  integral  in  equation  (8)  can  be  evaluated 
as : 


5 


i 


(x'-£')  +  r'  +  l-2r'cos 


CO 


(9) 


Thus , 


i_  PTw  _ x _  i 

p  ~  L  r,  .2  .2„  ,  ...  u 


[l  +  r*  +  x'-2r'  cos  9 1 


ac 


(10) 


It  should  be  noted  that  the  perturbation  pressure 
at  the  point  P(x,y,z)  or,  using  polar  coordinates, 
P(  x,r,a  ),  given  by  equation  (10),  is  that  due  to 
one  blade  at  azimuth  position  9  .  This  pressure 
is,  of  course,  dependent  on  the  angle  9'-9-a, 
which  represents  the  relative  position  between 
the  blade  and  the  point  under  consideration.  For 
a  given  point,  the  pressure  will  vary  as  the 
rotor  revolves  and,  hence,  is  a  time-varying 
quantity. 

The  time  dependency  of  the  pressure  at  any  point 
may  be  removed  either  by  using  the  concept  of  an 
infinite  number  of  blades,  or  by  averaging  the 
pressure  during  one  revolution  of  the  blade. 

This  is  equivalent  to  stating  that 


27 r 


27T 


(11) 


6 


Since  the  flow  induced  solely  by  helical  vortex 
is  irrotational  within  the  helical  cylinder, 
there  follows  that 


27 T 

J*r'ua^d#'=0  (12) 

o 

Also,  uqc  is  independent  of  9'  ,  so  that  the  last 
term  in  equation  (12)  is  equal  to  />curuac  •  The 
normal  average  induced  velocity  may  be  defined 
by 


27T 


un=27  /u Nd0' 


(13) 


Equation  (11)  therefore  becomes 


Uu  = 


N  27T/3U 


27 r 

J  p'dfi'-pwr'u^ 


(14) 


When  equation  (10)  is  substituted  into  equation 
(14), 


277- 

UN  =4ttU  I  2ir  J  T 


'ad' 


l  27r  0  [uAx'VcosS'] 


1/2 


(15) 
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t 


r\>|XD 


The  integral  in  equation  (15)  can  be  transformed 
into  the  complete  elliptical  integral  of  the 
first  kind  by  the  following  substitutions.  Let 


(3  = e ' 


2  2 


b  =(r'  +  I)  +x 
c2  =  (r'-l)2  +  x'  J 


(16) 


There  follows  that 


1 1  /+«'*-  2r'  cos  8'  -  -cg  1  cos  /3 

2  2 


.2  .2  (3  2  2  £  /  -i  7\ 

=  b  sin  —  +  c  cos  (17; 


Letting 
-  JL_S 

‘7 


x'd0' 


7T 

2 


=  4x'  f - 

i.fuZ 


dS 


27T 

^ ^l  +  r'  +x'  -2r'cos  0'j  cn/b^cos  2S  +  c2  sin  2  8 


7T 


_4_x 

b 


£ 

;/ 


dS 


(18) 


0  /.-k2  sin  2  8 
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i’here 


k 


2 


Substituting  equation  (18)  into  equation  (15) 
yields 


UN 


Tu> 
4  7T  U 


2% 

U^K(k) 
t r  b 


(19) 


where  K(k)  represents  the  complete  elliptical 
integral  of  the  first  kind.  For  N  blades  the 
normal  average  induced  velocity  is  obtained  by 
multiplying  the  right-hand  side  of  equation  (19) 
by  N. 

The  total  thrust  of  the  rotor  is  given  by 

R 

T  =  N  JpujrTdr^p^-rU  (20) 

o 

There  follows  that 


N  r  oi  _ _ _ 

47rU  2wR2  p\J 


(21) 


Substituting  equation  (21)  into  equation  (19), 
we  find  that 
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T 


(22) 


UN  s 


2  pir 


r  2x  1 

-5-  l+^-K(k) 

R  U  L  ^ 


When  equation  (22)  is  expressed  in  terms  of  uQ, 
given  in  equation  (2),  there  finally  results 


u 


N 


2x' 

1+  ^-K(k 
it  b 


) 


(23) 


(ii)  Tangential  Velocity  Component,  ua 

In  axial  flow  condition,  the  central  vortex  is 
shed  down  along  the  rotor  axis,  and  hence  the 
velocity  induced  by  the  central  vortex  is  always 
in  the  tangential  direction.  Referring  to 
Figure  3  and  utilizing  the  Biot-Savart  law,  we 
find  that  the  velocity , duQc  ,  induced  at  point 
P(  x,r,a  )  by  the  central  vortex  element  ds  can 
be  expressed  as  follows: 


T  dsxa 
du„  =T—  T= 


(24) 


x 

4r 


{ 


} 


e 


a 


where  ea  is  the  unit  vector  in  the  tangential 
direction  opposite  to  ou  . 
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By  integrating  equation  (24)  from  C  =  0  to 
£  =  00  ,  the  velocity  induced  by  the  central 
vortex  is  obtained  as 


4 irr  R 


,  ,2  ,2. .1/2 

(r  +x  ; 


+  1 


(25) 


For  an  N-bladed  rotor,  this  induced  velocity 
becomes 


-  Nr 

[c~  47rr'R 


,  ,2  ,2  1/2 

(r  +  x  ) 


(26) 


It  should  be  noted  in  equation  (26)  that  the 
velocity  induced  by  the  central  vortex, Uac  > 
is  independent  of  the  angular  position  of  the 
blade  and  the  point  P.  Its  time-averaged  value, 
UaC|  ,  can  therefore  be  represented  by  the  same 
expression. 

In  considering  the  total  time-averaged  tangential 
velocity,  ua  ,  the  effects  of  the  bound  vortex 
and  the  helical  vortex  can  be  neglected.  Ua  is 
therefore  given  by  the  same  expression  as  Uq 
Thus ,  c 


“a 


NT 

47rRr' 


x 


2  2  1/2 
(r‘  +  x*  ) 


+  1 


(27) 


(iii)  Radial  Velocity  Component,  Ur 

When  the  normal  and  tangential  velocity 
components  are  known,  the  radial  velocity 
component  is  derived  from  the  equation  of 
continuity,  which  states 


11 


i 


(28) 


t^-7 


<3ua 

dff 


<3un 

dx 


=  0 


The  minus  signs  indicate  that  ur  and  ua  are  in 
the  opposite  direction  of  r  and  6'  ,  respec¬ 
tively. 

As  noted  in  Reference  4,  the  normal  velocity 
component  uN  does  not  appreciably  vary  with 
radial  distance  r  ;  hence,  equation  (23)  can  be 
simplified  as  follows: 


uN  (r,,x,)=uN(',=0,x,)=u0[l  +  *  1 

L  Yu1** 


(29) 


When  equations  (27)  and  (29)  are  substituted 
into  (28)  and  with  the  boundary  condition 
ur  =  0  at  r  =  0,  we  find  that  the  time- averaged, 
radial-induced  velocity  is 


(30) 


b.  Instantaneous  Induced  Velocity 

( i )  Normal  Velocity  Component,  un 

The  instantaneous  value  of  the  normal  induced 
velocity  component  can  be  obtained  by  using  a 
method  similar  to  that  employed  for  the  average 
value.  However,  in  this  case  the  effect  of  the 
bound  vortex  must  be  retained. 
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Again  referring  to  equation  (3)  and  using 
equations  (6)  and  (7),  we  find  that  the 
perturbation  pressure  can  be  written  as 


P  =  "  4tt 


.2  .2  ,  1/2 
(1+  x  +r  -2 r  cos  9  ) 


-/>U(uNb)7>a>r'(uac4uab) 


The  total  normal  induced  velocity  is 


un  =  uNh+  uNb 


(32) 


Combining  equations  (3),  (31),  and  (32),  there 
follows  that 


CJ 


4ttU 


1  + 


(Ux'2+  r,2-2r'cos0')l/2 


(33) 


Due  to  the  fact  that  the  total  tangential 
velocity  component  is  small  compared  to  the 
normal  component  induced  in  the  rotor  wake,  and 
since  the  tangential  velocity  induced  by  the 
helical  vortex,  t  is  only  a  small  fraction 

of  the  total  tangential  velocity  induced,  u<*h 
in  equation  (33)  is  therefore  neglected  in  this 
analysis,  and  the  normal  induced  velocity 
becomes 


un  = 


feu 

4ttU 


1  + 


x 


(l+x,2+r,2-2r,cos0,)l/2 


(34) 
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In  the  above  equation,  the  normal  velocity 
component  induced  by  the  blade  bound  vortex, 
uN^  ,  can  be  derived  directly  from  the  Biot- 
Savart  law  as  follows: 

Referring  to  Figure  4,  we  find  that  the  velocity 
dub  induced  at  point  P  by  the  bound  vortex 
element  ds  is  given  by 


T  dsxa 

dUb  =^IiF 


(35) 


When  the  above  vector  quantities  are  resolved 
into  components  along  X,  Y,  Z  axes,  there  results 


diTb  =duxb'’+duyb"+duZb*  (36) 

dsxa=Rz(r'  sin  fi'ds'i-x'  sin  fids'  J-x'  cos  fids'k)  (37) 

lal3  =  R3  (x,%r,%*,2-2r,»,C05  9,)s/*  (38) 


where 


ds'  = 


and  i,  J,  and  k  are  unit  vectors. 


The  normal  velocity  component  can  be  obtained 
from  equations  (36),  (37),  and  (38)  by  solving 
for  the  i^  component  of  equation  (35);  thus, 
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du”b  =  47 T  R 


(39) 


_ r*  sin  ff'ds' _ 

. 2  .2  .2  .  .  3/2 

(x  +  r  +  s  - 2r s  cos  0  ) 

When  equation  (39)  is  integrated  from  blade  root 
to  tip,  there  results 

_ r'  sin  fl'ds1 _ 

/  |2  (2  ,  qi>3/2 

(x  +  r  +s  -2r  s  cos  8  ) 


UNb=  4ttR 


i 

/ 


Tr'  sin  9' 


-r'  cos  0' 


47rR(x'%r'2sin20')L(«2*r'2+:-2r'cos  e')'/2 


r'cos 
.27172 


'  .2 

(x  +r 


) 


(40) 


When  equation  (40)  is  substituted  into  (34),  the 
normal  velocity  induced  by  one  blade  in  the  rotor 
wake  is  obtained  as 


.  Tcu 


Un=4ttU 


I  +  2  2  1/2 

L  (Ux  +r'  -2r'cos  Q')  . 


T r'  sin  9' 

47rR(x'  ♦r'  sin  9') 


-r'cos  &  r'cos#1 

1/2  +  .2  .2  1/2 


[^l  +  x' +r'  -2r'cos#V~  (x‘*+r'~) 


(41) 


If  we  neglect  the  mutual  interaction  between  the 
vortex  systems  representing  different  blades,  the 
normal  velocity  induced  by  an  N-bladed  rotor  is 
obtained  simply  by  superposition  of  the 
velocities  induced  by  each  blade.  This  procedure 
is  expressed  mathematically  as  follows: 

for  an  N-bladed  rotor, 
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\ 


rVsinS'n 


4ttR(x,2+  r'2  sin  2  #'n  ) 


M-r'cos  0'n  r'cos  0'n\ 

- + - 

A  B 


(42) 


where 


A  =  (1  +  x'  +  r1  -  2r'cos  0'n  )  / 


,2  .2  1/2 

B  =  (x  +  r  ) 


(43) 

(44) 


(ii)  Tangential  Velocity  Component ,  Uq 

The  tangential  velocity,  ua  ,  induced  by  the 
rotor  is  obtained  by  considering  the  effects  of 
the  central  vortex  and  the  bound  vortex.  The 
effect  of  the  helical  vortex  on  the  induced 
velocity  is  small  in  the  tangential  direction 
and  hence  is  neglected.  The  velocity  induced 
by  the  central  vortex,  uQc  ,  is  readily  derived 
in  equation  (26).  The  tangential  velocity 
induced  by  the  bound  vortex  can  be  obtained  by 
an  analysis  similar  to  the  one  of  the  previous 
section.  Again  referring  to  Figure  4  and  the 
Biot-Savart  law,  the  Y  and  Z  components  of  the 
velocity  induced  by  the  bound  vortex  are 
obtained  as 


F  x'  sin  9 
4  ttR 


f  ds' 

J  ,  ,2  >2  >2  0  ■  -  ni,3/2 

o  (x  +r  +  s  -2r  s  cos  6  ) 


(45) 


If 


r  x'cos  6 
Jzb  4ttR 


r  ds 

J  ,  ,2  ,2  ,2  ,  , 

n  (x  +  r  +  s  -2rs 


cos  @')3/2 


(46) 


Resolving  Uy^  and  uz^  in  the  tangential 
direction,  there  results 


Ufltb  =uyb  sin  a  +  Uzb  cos  a 


(47) 


When  equations  (45)  and  (46)  are  substituted 
into  (47)  and  integrating,  the  tangential 
velocity  induced  by  the  bound  vortex  on  one 
blade  becomes 


rx'cos#'  f  l-r'cos  0' .  r'cos#' 


“b  4wR  (x'Z+  r,Z  sin2#1  ) 


_  LL££ 

i1  \ .  A 


B 


(48) 


For  an  N-bladed  rotor, 


N 

u<*b=£ 

n=l 


Tx'costf'n  /l-r'cos  &'n  r'cosfl'n 

L4.R(x'2*rZs,nVn)V  A  *  ~B~ 


(49) 


When  equations  (26)  and  (49)  are  combined  the 
total  tangential  velocity  component  becomes 


NT  /  x1  \  £ 
u<2= - 1  -  +  l)  ~L 

4-rrRr  \B  /  n=l 


rx'cosg'n  ( l-r'cos  g'n 

2  2  2  1  a 

4ttR(x'  +r'  sin  0'n)\ 


r'  cos0'n 

B 
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(50) 


(iii)  Radial  Velocity  Component, 

From  equations  (45)  and  (46)  it  can  be  noted 
that  the  radial  velocity  induced  by  the  bound 
vortex  is  zero,  since 


u»'b  =  uZbSina~uyb 


cos  a  =  0 


Since  the  velocity  induced  by  the  central 
vortex  is  only  in  the  tangential  direction,  the 
radial  velocity  component  is  induced  only  by  the 
helical  vortex.  Assuming  furthermore  that  uN 
does  not  vary  appreciably  with  r1  ,  equation  (42) 
reduces  to  the  averaged  value  given  in  equation 
(29).  There  follows  that  the  radial  induced 
velocity,  ur  ,  is  the  same  as  its  time-averaged 
value  given  in  equation  (30). 


ur  =  ur 


up  f  r1 

2  Lx'*)3* 


(51) 


2.  Triangular  Disk  Loading 

Based  on  the  results  of  uniform  loading,  the  normal  induced 
velocity  for  a  triangularly  loaded  rotor  is  obtained  by  super 
position  of  uniformly  distributed  loadings.  In  this  case, 
the  circulation  of  the  bound  vortex  along  the  rotor  blade, 
instead  of  being  constant,  will  vary  along  the  length  of  the 
blade,  and  helical  vortices  will  be  shed  all  along  the  blade, 
thus  forming  a  vortex  sheet  of  helicoid  shape.  This  vortex 
sheet  is  generated  in  a  way  similar  to  the  flat  vortex  sheet 
behind  a  lifting  wing  of  nonuniform  loading.  In  a  tri¬ 
angularly  loaded  rotor,  this  helical  vortex  sheet  replaces 
the  blade  tip  helical  vortex  and  the  central  vortex  which 
exist  in  the  case  of  uniform  loading,  as  previously 
discussed. 

a.  Time-Averaged,  Induced  Velocity 
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(i)  Normal  Induced  Velocity  Component,  un 

This  velocity  component  of  a  triangularly  loaded 
rotor  is  obtained  from  the  results  derived  in 
uniform  loading  by  a  method  of  superposition. 
Figure  5  is  a  graphical  representation  of  a 
triangularly  loaded  blade  with  circulation, 
r=rt(r/R)  .  The  blade  tip  circulation,  Tf  ,  may 
be  evaluated  since 


T=  j  Npu)Tdr=  N^>a»rf^- 
o 


Consequently, 


rt 


3T 

N/owR2 


(52) 


It  can  be  seen  in  Figure  5  that  the  velocity 
induced  by  the  triangular  loading, r=rf(r/R)  , 
can  be  obtained  by  the  difference  between  that 
due  to  the  assumed  uniform  loading,  Tf  ,  shown 
in  dotted  line  in  Figure  5,  and  that  due  to  the 
integrated  loading  of  dTi  from  r  =  0  to  r  =R  . 

The  normal  average  induced  velocity  of  an 
N-bladed  rotor  induccsd  by  the  uniform  loading, 
Tf  ,  it  point  P(x,r)  is  obtained  from  equation 
(29). 


t 


n  r  _  (u 

4  ttU 


(53) 


Similarly,  the  normal  velocity  induced  by  the 
differential  uniform  loading,  d T,  ,  at  point  P 
can  be  expressed  as  follows: 
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i 


Substituting  equations  (55)  and  (56)  into 
equation  (54),  we  have 


It  may  be  seen  (Reference  6)  that  the  velocity 
Induced  by  a  vortex  cylinder  is  discontinuous 
across  the  vortex  boundary  and  that  outside  the 
boundary  the  induced  velocity  is  negligible. 
Thus,  the  velocity  induced  at  a  given  radial 
station  is  affected  only  by  vortices  shed  by  the 
blade  outside  this  station.  The  induced 
velocity,  u^j  ,  at  point  P(x,r)  due  to  the 
uniformly  distributed  loadings,  dTj  ,  super¬ 
imposed  along  the  blade  is  therefore  obtained  by 
integrating  equation  (57)  from  Rj=r  to  Rj  =R  . 
Thus , 
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The  normal  velocity,  uN  ,  induced  by  the 
triangular  loading,  T  ,  is  therefore  obtained 
from  equations  (53)  and  (59)  by  the  relation 
UN=uNt"uNj-  Thus, 


_  nTjcu 
U"=4^u 


■N 

► 


(60) 


Substituting  equation  (52)  into  equation  (60), 
there  results  that 


3T 


/,  , 


(61) 


Expressing  uN  in  terms  of  u0  ,  equation  (61) 
becomes 
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(62) 


qn  =  Ygo^ 


l  I 

r  +  x 


•  •v/i  +  x' 2 


+  \/U 


.2  -1 


,  /r^2 

r  +v  r  +  x 


(ii)  Tangential  Velocity  Component, 

In  a  triangularly  loaded  rotor,  the  helical 
vortex  sheet  replaces  the  central  vortex  which 
exists  in  the  case  of  uniform  loading.  In 
obtaining  the  tangential  velocity  component  in  a 
triangularly  loaded  rotor,  the  approach  utilizing 
the  effect  of  a  central  vortex,  considered  in  the 
case  of  uniform  loading,  is  therefore  not 
applicable.  This  velocity  component,  however, 
can  be  obtained  by  numerically  averaging  the 
instantaneous  tangential  velocity  induced  by  the 
bound  vortex.  This  instantaneous  value  is  given 
by  equation  (72). 

(iii)  Radial  Velocity  Component, Ur 

This  velocity  component  is  derived  similarly  to 
that  in  the  case  of  a  uniformly  loaded  rotor. 
Since  the  average  tangentially  induced  velocity, 
ua  ,  is  independent  of  9',  the  radial  velocity 
component,  ur  ,  is  obtained  by  substituting  from 
equation  (62)  into  the  equation  of  continuity, 
equation  (28).  After  integrating  and  with  the 
substitution  of  the  boundary  condition,  ur  =  0 
at  r  =  0,  this  velocity  is  obtained  as 


u 


r 


+  1 


vA-(l  +  v4*  x'  ) .  + 


(63) 
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b.  Instantaneous  Induced  Velocity 

(i)  Normal  Velocity  Component ,  un 

The  instantaneous  value  of  this  velocity 
component  is  obtained  by  the  same  superposition 
method  utilized  in  deriving  the  time-averaged 
value.  Referring  again  to  Figure  5,  the 
instantaneous  normal  velocity, uNj  ,  is  obtained 
from  equation  (42)  as 


Tt  r'sin0'n 
47tR(x'  +r'Zsin20'n) 


rl-r'cos0'n 

r'cosS'nlN 

L  A 

+  B  Jj 

(64) 


The  normal  velocity  induced  by  the  incremental 
uniform  loading,  dTj  ,  is  similarly  obtained  as 


u 


Nj 


Tt  rWnd(^) 

y.  D/  |2  ,  2—| 

47rR{x  +r  sin  on) 
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It  should  be  noted  in  the  above  equation  that 
the  second  part  is  the  velocity  induced  by  the 
blade  bound  vortex.  This  velocity  is  not 
affected  by  the  vortex  boundary;  hence,  the 
integration  is  performed  from  blade  root  to 
tip  where  the  limits  are  Rj  =  0  and  Rj  =  1. 

R  R 

After  integrating  equation  (65),  there  results 


u 


H  r'sinfl'n 

+ - 

47rR  (  x'  ♦  r'  sin2#^) 


(66) 


where 


l-r'cos0'n  +  A 

r'  ^  l-cos@'nj+  x  +2r'  ^l-cos  0'nj 


(67) 
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The  resultant  instantaneous  normal  induced 
velocity, uN  ,  is  therefore  obtained  from 
equations  (64)  and  (66)  by  the  relation 
uN=uNt-uNj-  Thus, 


(ii)  Tangential  Velocity  Component ,ua 

In  deriving  the  instantaneous  value  of  the 
tangential  induced  velocity,  it  will  be  assumed 
that  the  effect  of  the  helical  vortex  can  be 
neglected.  Furthermore,  in  the  case  of  a 
triangularly  loaded  rotor,  there  exists  no 
central  vortex.  Hence,  the  instantaneous 
tangential  velocity  component,  ua  ,  is 
obtained  by  considering  only  the  effect  of  the 
bound  vortex.  Referring  to  Figure  5  and 
equation  (49),  the  tangential  velocity  induced 
by  the  uniform  loading, Tt  ,  at  point  P(x,r  ) 
can  be  obtained  as  follows: 


Tf  x'cos0'n  rl-r'cos@n 


47tR(x'  +  r'  sin20'n) 


-i 

il  \ 


(69) 


25 


The  tangential  velocity  induced  by  the  incre 
mental  uniform  loading,  dTj  ,  is  similarly 
obtained  to  be 


n  j.  i  r,  «Wrfl(^)r  (  -r'cosS'n 
V  '^rRi 


e'nl  ([r]  «,2*r'2-2f,(|)cose,l „) 


1/2 


r'cosfl'n 
(x,2+  r,Z  )  1/2 


]} 


(70) 


The  integration  of  equation  (70)  yields 


v-i[ 

n-l  L 


rtx'cos0'n  /r'cos9'n 


n=lL47rR  (x'2+r'2sin2^) 


/  rcosbn  \  1 

(71) 


The  resultant  instantaneous  tangential 
velocity,  ua  ,  is  therefore  obtained  from 
equations  (69)  and  (71)  from  the  relation 


ua  =  uat  ~uoj  •  Thus- 


~s 


rtx'cos0'n  ^l-r'cos0'n 


X  l  ,2  I2  2rtl  \ 

n=l  4?rR^x  +r  sin  9n ) 


-A+  B 


)] 


(72) 


(iii)  Radial  Velocity  Component,  ur 

The  method  of  superposition  of  uniform 
loadings  is  again  utilized  in  deriving  this 
induced  velocity  component  of  a  triangularly 
loaded  rotor.  The  radial  velocities  induced 
by  the  uniform  loading, Tt  ,  anddFj  ,  at  point 
P(x, r  )  can  be  obtained  from  equation  (51) 
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as  follows: 


The  integration  of  equation  (74)  yields 


I 


Subtracting  equation  (75)  from  (73),  the 
resultant  instantaneous  induced  radial 
velocity  component ,  ur  ,  becomes 


B.  YAWED  FLOW 

The  rotor  in  yawed  flow  is  shown  in  Figure  6.  It  is  assumed 
that  the  vortices  shed  from  the  rotor  in  the  yawed  flow 
condition  are  shed  with  a  velocity  equal  to  the  resultant 
averaged  velocity  of  the  flow  through  the  rotor.  This 
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resultant  velocity  may  be  represented  by  a  component, 

Vcos<£+u,  perpendicular  to  the  rotor  plane  and  the 
component ,  Vsin  <p0  ,  parallel  to  the  plane.  As  shown  in  Figure 
6,  the  rotor  wake  is  therefore  swept  back  at  an  angle  X  with 
respect  to  the  rotor  axis,  the  value  of  X  being 


Vsin</> 

v  =  tan  1  - 7 — — 

A  Vcos  </>  +  u0 


(77) 


The  resultant  normal  velocity  at  the  rotor  plane,  U  ,  in 
yawed  flow  condition  is  given  by 


U  =  V  cos<£+  u0 


(78) 


where 


Zp  A^/(  Vcos^>  +  u0)2  +  (Vsin<£) 


(79) 


For  any  given  condition  of  flight,  V  and  <£  are  specified; 
hence,  u0  and  U  can  be  obtained  by  the  simultaneous  solution 
of  equations  (78)  and  (79). 

The  coordinate  system  defining  the  induced  velocities  in 
yawed  flow  condition  is  shown  in  Figure  7.  Xi,  Yi,  and  Zl 
are  the  axes  of  a  coordinate  system  with  the  X^-Z^  plane 
rotated  an  angle  X  with  the  X-Z  plane  such  that  the  X-^-axis 
coincides  with  the  axis  of  the  swept-back  wake.  The  angular 
position  of  point  P  is  defined  by  any  of  the  following  three 
angles,  respectively: 

(i)  a  =  angle  defining  point  P  with  vertex  on 
X-axis,  measured  from  positive  Y-axis 
parallel  to  the  Y-Z  plane. 


28 


(ii)  a  =  angle  defining  point  P  with  vertex  on 
Xi-axis,  measured  from  positive  Y-axis 
parallel  to  the  Y-Z  plane. 

(lii)  a,  =  angle  defining  point  P  with  vertex  on 
X^-axis  parallel  to  the  Yi-Z^  plane. 

All  angles  are  positive  in  the  direction  of  .  The  trans¬ 
formation  equations  between  the  X]_,  Yi,  Zi  and  X,  Y,  Z 
systems  are  given  as  follows: 


x(  =  x -z  tan  X 

y  =  -1- 

y'  cos  X 

z,  =  Xtan  X  +  z 

rM^xf  +  («'onX”)2 

tan  a (  =  sinX+ cosxj  tana 

f  z  +  x  tan  X  1 
tan  a  =  -  - y - J 


(80) 


The  normal  velocity  component  induced  by  a  lifting  rotor  in 
the  yawed  flow  condition  is  derived  from  the  results  of  the 
equations  in  axial  flow  with  an  additional  parameter  X 
defined  by  equation  (77).  The  vortex  system  assumed  in  this 
case  is  essentially  the  same  as  that  in  axial  flow,  but 
located  in  a  differenct  position  with  respect  to  the  rotor. 
The  central  vortex  is  now  shed  along  the  Xi-axis,  and  the 
helical  vortex  is  shed  down  from  the  blade  and  winds  around 
the  X^-axis.  The  bound  vortex  is  fixed  to  the  blade  and 
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hence  remains  in  the  same  position  as  in  the  case  cf  the 
axial  flow  condition.  The  velocities  induced  by  these  three 
line  vortices  are  discussed  separately  as  follows: 

a.  Normal  Velocity  Induced  by  the  Trailing  Helical 
Vortex 

It  has  been  shown  in  Reference  4  that  the  normal 
velocity  induced  by  the  trailing  vortex  at  point  P, 
defined  by  the  coordinates  ( x',  ,  ,  a  )  in  the 

slipstream  of  a  rotor  yawed  at  an  angle  %  ,  is 
equivalent  to  tho  induced  velocity  existing  at  the 
corresponding  point  (  x' =  x',  ,  r '  =  r',  ,  ara  ) 

in  an  unyawed  rotor  multiplied  by  cos  X .  This  can 
be  expressed  as 


Wi  ,5)x=x=cosX[Uxh(x=x',r=ri  ,a=a)x-o 


where 


R2  =  R2  ( I  -  sin2  a  sin2  X  ) 


(81) 


(82) 


(83) 


According  to  equation  (81)  and  given  the  coordinates 
X|  ,  r | ,  and  a  of  a  given  point  P  in  the  rotor  wake 
and  the  wake  angle  X ,  the  normal  velocity  induced 
by  the  helical  vortex  of  a  yawed  rotor  at  that  point 
can  therefore  be  obtained  fron  the  expression  derived 
for  the  axial  flow  condition. 
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b .  Velocity  Induced  by  the  Central  Vortex 


In  axial  flov.’  condition,  the  central  vortex  is  shed 
along  the  X-axis.  The  velocity  induced  by  the 
central  vortex  is  therefore  in  the  tangential 
direction  and  does  n  t  contribute  to  the  normal 
velocity  induced  in  the  rotor  wake.  In  yawed  flow 
condition,  however,  th-  central  vortex  is  shed  along 
the  Xi-axis,  and  the  vt  jcity  induced  has  a  component 
along  the  normal  direction  or  the  X-axis.  This 
component  is  given  as  follows: 


uw  =u„  .cosa.sinX  (84) 

™c  ai  1 


where  ua,  is  the  total  velocity  induced  by  the 
central  vortex,  and  uN  is  the  velocity  component 
along  the  X-axis. 

The  total  induced  velocity,  UQ|  ,  at  point 

P(  r'j.x',  ,  a  )  can  be  expressed  similar  to  equation 

(26),  as  follows: 


n  r*  r 


47rRr, 


l'2  '2\ 

I*.  * r,  j 


1/2 


+  I 


(85) 


Substituting  equation  (85)  into  (84)  yields  the 
normal  component  of  uQ|  along  the  X-axis: 


n  T cos  a,  sin  X 
4  7r  R  r! 


(fo 


(86) 
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where 


c .  Velocity  Induced  by  the  Bound  Vortex 

Since  the  position  of  the  bound  vortex  is  not 
affected  by  the  yawed  flow,  the  bound  vortex-induced 
velocity  derived  in  axial  flow  condition  is  directly 
applicable  in  yawed  flow  condition. 

From  the  above  consideration,  the  normal  velocity  ir  ,  />d 

at  any  point  in  the  wake  of  a  rotor  in  yawed  flight  condition 
can  be  obtained  from  those  expressions  derived  for  the  axial 
flow  condition  by  introducing  a  new  set  of  position  parameters 
defining  that  point  and  the  yawed  wake.  These  parameters  are 
X|  ,  r',  ,  a,  ,  a  ,  and  X  .  For  a  given  condition  of  flight, 
these  parameters  can  be  readily  calculated  by  equations  (77), 
(80),  and  (82). 

1 .  Uniform  Disk  Loading 

a.  Time -Aver aged.  Normal  Induced  Velocity  un 

According  to  equations  (81)  and  (22),  the  time- 
averaged,  normal  induced  velocity,  uN ,  can  be 
expressed  as  follows: 


u 


N 


Tcos  X 
2/jttR2U 


(87) 


where  bi  and  kp  are  given  in  equations  (16)  and  (18) 
in  terms  of  x',  and  r . 

From  equations  (77)  and  (78),  cos  X  can  be  expressed 
as 


32 


u 


(88) 


cosX 


v^(vsin<^)  +  (vcos<£+ 


Substituting  equation  (88)  into  (87),  we  have 


2/3  7rR2v/(Vsin<^}2+(Vcos<^  ♦  uQ)2 


l4^r  K<k!> 

L  7Tb.  1  J 


(89) 


Substituting  u0  Qf  equation  (79)  into  equation  (89) 
and  utilizing  the  approximation,  uN  (  r]  ,  x‘,  )  = 

uN  (  x\  -  0,  x',  ),  equation  (89)  becomes 


u 


N  = 


n  M,/2 

(1+  *,) 


(90) 


b.  Instantaneous  Normal  Induced  Velocity 

In  the  axial  flow  condition,  this  velocity  component 
is  given  by  equation  (42).  In  this  equation,  the 
first  part  gives  the  velocity  induced  by  the  helical 
vortex,  whereas  the  second  part  gives  the  velocity 
induced  by  the  bound  vortex.  In  the  yawed  flow 
condition,  the  bound  vortex  effect,  uN^  ,  remains 
the  same.  Thus, 


N 

UNb=I 

— 


f  r'sin@'n 


"='  4wR(*'a*r'25inV) 


l-r'cos0'n  r'cos0'n 
A  +  B 


(91) 
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The  normal  velocity  induced  by  the  helical  vortex, 
,  however,  becomes: 


N 


n=l 


rcucosX  ’  _Xj_ ' 

4ttU  L  A,  . 


(92) 


where 


A.  =  (I  +  x'2+  r'2-2r  cos#n) 


n 


and 


®n  'V5 

In  addition,  the  normal  component  of  the  central 
vortex- Induced  velocity,  ukc  ,  given  in  equation  (86) 
must  be  included.  The  total  normal  induced  velocity 
in  yawed  flow  condition,  uN  ,  can  therefore  be 
obtained  as  follows: 


''n=X  +  1Vu»c 


n  TwcosX  r  x'|  ]  r r'sin0'n  r  l-r'co$0'n 
un  Z  /iwM  I1  aJ  47rR^x2+r2sin20'nH  A 


ri=l  4^U  L 


r'cos0'nl  NfcosajSinX  r  x', 


B 


1  Ni  cosa,  sin  a  r  x,  ~\ 

- 3— 5“i —  F  +  l  (93> 

J  47rRr,  LB.  J 
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2.  Triangular  Disk  Loading 


The  velocity  induced  by  a  triangularly  loaded  rotor  in  yawed 
flight  condition  can  be  obtained  from  those  results  derived 
for  axial  flight  condition  in  exactly  the  same  way  as  in  the 
case  of  the  uniformly  loaded  rotor  previously  discussed. 
However,  the  normal  velocity  induced  by  the  central  vortex, 
unc ,  is  neglected  in  this  case  due  to  the  fact  that  for  a 
triangularly  loaded  rotor  there  no  longer  exists  a  central 
vortex.  The  time-averaged  and  instantaneous  normal  induced 
velocity  components  are  thus  obtained  from  equations  (62) 
and  (68),  respectively,  and  are  presented  below. 

a.  Time -Averaged,  Normal  Induced  Velocity,  un 


where  uQ  is  given  in  equation  (79). 
b.  Instantaneous  Normal  Induced  Velocity,  un 


n=l  ' 


rtr'sin#'n 


4ttR 


(*' 


2  ,2  2 

+  r  sm 


•  \  - 


I-  r'cos0'n 


®r> ) 


-  A  +  B 


(95) 
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where 


l-r',cos0n  +  A, 

r'i  ( l-cos0n)  ♦  (  x'f  +  2rf[l-cos0n]  )'/2 
C.  FUSELAGE  EFFECTS 


(96) 


The  effect  of  the  fuselage  on  the  induced  velocities  inside 
the  rotor  wake  is  derived  here  for  uniformly  loaded  and 
triangularly  loaded  rotors  in  the  axial  flow  condition.  The 
fuselage  is  represented  by  s  three-dimensional  doublet  placed 
in  the  nonuniform  flow  field  corresponding  to  the  time- 
averaged,  axial-induced  velocity  in  the  rotor  wake.  Using 
this  flow  pattern,  the  strength  of  the  doublet  is  determined, 
and  the  zero  streamline  of  the  flow  corresponding  to  the 
fuselage  shape  is  established.  The  velocity  at  a  given  point 
in  the  rotor  wake  can  be  computed  by  the  superposition  of  the 
velocity  induced  by  the  rotor  and  that  induced  by  the  doublet 
of  strength  M.  In  the  case  of  a  triangularly  loaded  rotor, 
there  is  an  upwash  velocity  region  near  the  rotor  axis  and 
immediately  below  the  rotor  plane,  as  predicted  by  the  flow 
field  analysis  of  a  free  rotor.  The  present  analysis  of 
the  fuselage  effect  is  applicable  only  when  the  fuselage  is 
placed  within  this  upwash  region. 

1,  Uniform  Loading 

The  t ime- averaged ,  normal  induced  velocity  in  the  wake  of  a 
free  rotor  is  given  in  equation  (29)  as  follows: 


The  flow  represented  by  equation  (97)  in  the  rotor  wake  is 
axisymmetric  with  respect  to  the  rotor  axis.  The  stream 
function  for  such  a  flow  is  defined  as 
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(98) 


Substituting  equation  (97)  into  (98)  yields 


u 


o 


I 


<wn- 


2-rrR  r  \  dr' 


(99) 


Transforming  equation  (99)  into  the  coordinate  system 
fixed  to  the  fuselage  center  of  gravity  shown  in  Figure  8 
yields 


(100) 


where 


Integrating  the  above  equation,  we  have 


(101) 
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Equation  (101)  represents  the  stream  function  of  the  induced 
flow  inside  the  wake  of  a  free  lifting  rotor.  The  fuselage 
is  now  replaced  by  a  three-dimensional  doublet  located  at 
the  center  of  gravity  of  the  fuselage.  Its  stream  function 
is  given  by 


M  (  sin2/3\ 
2R  '  s’  / 


(102) 


where 


s 


_s_ 

R 


and  s  is  the  radial  distance  measured  from  the  fuselage 
center  of  gravity  to  the  point  in  space. 

Combining  equations  (101)  and  (102)  gives  the  total  stream 
function  inside  the  wake.  Thus, 


•  ,  ,  §  M  /  sin2/3\  *  .2  T . 

♦-V+0-2R  (— )+,rR  uor  L,+ 


(h'.x'l 


‘J\*  (h1  +  x1  )2 


(103) 


The  zero  streamline  of  the  flow  is  obtained  from  equation 
(103)  letting  ^  =  0,  which  yields 


-  ttR2G0 


(h'+  s'cos/3) 


(104) 
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It  is  found  in  the  above  equation  that  s  varies  only  slightly 
with /3  .  Equation  (104)  can  therefore  be  simplified  by 
letting  cos f3  =  0.  Thus, 


(105) 


For  an  assumed  spherical  fuselage  of  radius  a,  represented 
by  the  above  zero  streamline,  the  strength  M  of  the  three- 
dimensional  doublet  is  calculated  from  equation  (105)  by 
the  substitution  S’  =  a'  ,  where  a'  =  a/R. 


M  =  27rR3a’3u0( 


(106) 


Substituting  equation  (106)  into  (104)  yields  the  zero 
streamline  or  the  approximated  fuselage  shape  as  follows: 


cos/3  =—■ 
s 


(107) 


where 


K  =1  + 
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The  normal  velocity  induced  by  the  doublet  of  strength  M 
given  in  equation  (106)  can  be  expressed  in  terms  of  its 
stream  function  as  follows: 


I 

27rR2r'  dr' 


(108) 


Substituting  equation  (102)  into  (108)  yields 


(109) 


For  a  given  fuselage  of  equivalent  radius,  a,  the  doublet 
representing  this  assumed  fuselage  is  defined  by  its 
strength  M,  given  in  equation  (106).  The  normal  velocity 
induced  in  the  * otor  is  obtained  simply  by  adding  the 
velocity  induced  by  the  free  rotor  and  that  induced  by  the 
doublet  given  in  equation  (109).  Though  equation  (109)  is 
derived  from  time -averaged,  induced  velocity,  it  is  assumed 
to  be  applicable  also  for  instantaneous  velocity  induced 
inside  the  rotor  wake. 


Triangular  Loadini 


As  previously  stated,  this  analysis  of  the  fuselage  effect 
for  the  case  of  triangular  loading  is  applicable  only  when 
the  fuselage  is  placed  within  the  upwash  region  near  the 
rotor  axis.  This  is  due  to  the  fact  that  the  zero  streamline 
of  the  combined  flow,  thac  due  to  che  free  rotor  and  the 
assumed  doublet,  is  no  longer  a  closed  surface  where  reverse 
flow  exists.  The  fuselage  chape  therefore  cannot  be  estab¬ 
lished  in  this  case.  However,  if  the  equivalent  radius  of 
the  fuselage,  a,  is  within  the  upwash  region,  the  fuselage 
effect  can  be  analyzed  similar  to  the  case  of  uniform 
loading. 
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The  normal  average  velocity  induced  by  a  triangularly  loaded 
rotor  is  given  in  equation  (62)  as 


/  -j 

1+  V  I  ♦  * 


l  ,  fT~ 

r  +  v  r  + : 


(110) 


where 


x  =  (h'+  x  ) 


The  upwash  region  can  be  easily  determined  by  computing  a 
value  of  r'  from  equation  (110)  corresponding  to  the 
condition  ux  =  0  at  x‘  =  0.  If  a'<  the  computed  value  of 
r'  ,  then  the  f  blowing  analysis  is  applicable. 

The  combined  stream  function  of  the  flow  induced  by  the  free 
rotor  and  the  doublet  located  on  the  rotor  axis  at  the 
fuselage  center  of  gravity  is  derived  similar  to  that  of 
uniform  loading  and  is  obtained  as  follows: 


(111) 


41 


For  ^  =  0,  x  =  0,  and  r'  =  a'  ,  the  strength  of  the  doublet> 
M,  is  obtained  graphically  from  equation  (111).  Once  the 
strength  of  the  doublet  is  defined,  the  velocity  induced 
by  the  doublet  is  calculated  from  equation  (109),  and  the 
total  induced  velocity  at  a  given  point  in  the  rotor  wake  is 
computed  by  the  superposition  of  the  velocity  induced  by 
the  rotor  and  that  induced  by  the  doublet. 

D.  GROUND  EFFECT 

A  lift  device  operating  in  the  proximity  of  the  ground 
experiences  a  reduction  in  the  induced  velocity.  An  indica¬ 
tion  of  the  reduction  of  the  magnitude  of  the  time -averaged, 
normal  induced  velocity  is  shown  in  test  data  such  as  those 
in  References  7  and  8.  In  order  to  account  for  this 
reduction,  a  correction  factor,  kg,  has  been  obtained 
empirically  by  use  of  the  data  given  in  the  above  two 
references.  This  factor  is  obtained  as 


-  2  — 

kg  =1  -  0.9e  *  R 


(112) 


The  induced  velocity  in  ground  effect  can  then  be  computed 
from  the  out-of -ground  value  as  follows: 


UIGE=  kgU0GE 


(113) 


E.  WAKE  CONTRACTION 

The  analysis  presented  above  incorporates  the  assumptions  that 
the  vortices  shed  at  any  point  from  the  blade  remain  at  a 
constant  strength  with  time  and  form  a  helical  path  on  a 
cylinder  of  constant  radius.  The  latter  assumption  implies 
that  the  effect  of  wake  contraction  is  neglected.  Wake 
contraction  is  only  of  minor  importance  in  high  forward  or 
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climbing  speed  conditions,  but  it  is  a  significant  factor 
in  the  hovering  flight  condition. 

One  method  of  including  the  effect  of  rotor  contraction  for 
the  hovering  flight  condition  consists  simply  of  assuming 
that  the  induced  velocities  at  any  radial  location  r/R  can 
be  computed  on  the  basis  of  the  local  wake  radius.  Re,  i.e., 
u(r)corr  =  u(£  )*  A  more  l°g*-cal  method  of  accounting  for 

wake  contraction  is  illustrated  in  Figure  9.  The  induced 
velocities  at  any  point  experience  the  strongest  effect  by 
those  vortices  which  are  nearest  to  them.  Let  it  therefore 
be  assumed,  in  order  to  obtain  the  velocities  at  a  point 
Pl(x,/R, ,  r(/R),  that  the  outer  radius  of  the  cylinder  which 
incloses  the  shed  vortices  is  equal  to  that  of  the  local 
wake  radius,  Re(X,).  Referring  again  to  Figure  9,  this  then 
implies  that  the  distance  between  the  point  Pi  and  A,  located 

at  an  axial  position  smaller  than  Xt  ,  is  larger  than  that 
used  in  the  calculations.  On  the  other  hand,  the  distance 
between  the  point  Pi  and  the  vortex  element  B,  located  at  an 
axial  position  longer  than  X  (  ,  is  smaller  than  that  used  in 
the  calculations. 

The  induced  velocities  can  then  be  computed  using  the 
previous  theory  by  replacing  R  by  Re,  i.e., 


(  X| 

r_\\  / 

xi  ri  \ 

1  R  ' 

1  R  /CORRECTED- U\ 

Re’Re  ' 

(114) 


The  values  of  Re  can  be  obtained  from  experimental  data, 
such  as  those  presented  in  Reference  9. 
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III.  THE  TEST  PROGRAM 


Experimental  data  were  obtained  for  both  instantaneous  and 
time-averaged,  induced  velocit'  s  for  several  flight 
conditions  of  a  full-scale  helicopter.  The  details  of  the 
test  program  are  described  in  this  section. 

A.  TEST  HELICOPTER 

The  aircraft  utilized  for  the  test  program  was  a  commercial 
Hughes  helicopter,  Model.  269-A  (Army  designation  TH-55A). 

This  helicopter  was  chosen  for  this  program  for  the  following 
reasons : 


1.  Low  vibration  and  noise  levels. 

2.  Adequate  cabin  room  for  recording  instrumentation. 

3.  Unobstructed  360°  visibility  in  cabin. 

4.  Compact  in  size. 

5.  Acceptable  handling  qualities  in  hovering  and 
forward  flight. 

6.  Small  fuselage  size  to  minimize  fuselage  effects. 

Table  I  presents  the  principal  dimensions  and  physical 
characteristics  of  the  Hughes  269-A  test  helicopter. 

B.  TEST  INSTALLATION 

The  test  rig  utilized  to  measure  rotor  velocities  is  shown 
mounted  on  the  test  helicopter  in  Figure  10.  The  test  rig 
consists  of  a  rectangular  frame  with  its  vertical  members 
parallel  to  the  rotor  shaft.  Attached  to  the  vertical 
members  of  the  frame  are  two  horizontal  booms  (one  on  each 
side  of  the  helicopter) ,  upon  which  are  mounted  the  probes 
which  measure  time-averaged  and  instantaneous  induced 
velocities . 
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1.  Support  Structure 


The  boom  support  structure,  a  space  framework,  was  attached 
to  the  helicopter  at  existing  hard  points,  i.e.,  the  rotor 
mast  shelf,  skid  gear  cross  beam,  and  the  tail  cone  fuselage 
cluster,  with  a  minimum  amount  of  alteration  to  the  existing 
structure. 

The  rotor  mast  shelf  was  backed  up  with  an  additional  shelf 
structure  that  was  attached  by  picking  up  existing  rivet 
patterns.  This  shelf  backup  supported  the  upper  portion  of 
the  forward  vertical  columns.  The  lower  supports  for  these 
columns  were  provided  by  modifying  the  step  structure  which 
replaced  the  existing  steps.  The  aft  support  for  the  frame¬ 
work  was  accomplished  with  a  clamp  collar  that  fitted  the 
tail  boom  at  the  fuselage.  The  structure  installation  is 
shown  in  Figure  11. 

2.  Instrument  Support  Booms 

The  velocity  probe  booms  were  attached  to  the  vertical 
support  columns  with  clamp  fittings  to  allow  for  changes  in 
azimuth  and  height  positions.  The  booms  were  fixed  in 
azimuth  position  with  additional  tube  structure,  as  shown  in 
Figure  11.  The  booms  span  approximately  90  percent  of  the 
rotor  radius.  The  natural  frequency  of  the  booms  was 
measured  to  assure  that  it  did  not  coincide  with  the  rotor 
harmonics.  The  maximum  vertical  test  height  of  the  probe 
from  the  rotor  blades  was  based  on  a  blade  tip  to  tail  boom 
clearance  of  12.0  inches,  as  established  by  structural 
considerations.  The  horizontal  support  booms  could  be  moved 
aft  of  the  rotor  shaft  and  fixed  in  the  required  position 
for  recording  data  in  skewed  wakes  generated  at  forward  speeds. 

C.  INSTRUMENTATION 

A  block  diagram  of  the  instrumentation  system  utilized  to 
obtain  the  required  data  is  shown  in  Figure  12.  This  system 
includes  data  sensing,  recording,  and  playback  equipment. 

1 .  Data  Sensing  Equipment 

a.  Time-Varying  Velocities 
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Two-dimensional  nX"  configuration  quartz-coated 
cylindrical  hot  film  sensors  (mounted  on  6-inch- 
stem-length  probes)  were  utilized  to  sense  the  rotor 
wake  time-varying  induced  velocities.  The  principle 
of  operation  is  the  same  as  that  of  a  hot  wire 
sensor. 

The  hot  film  sensors  with  a  .006-inch  diameter  and  a 
.125-inch  length,  Model  TSI-DYN,  were  made  by  Thermo- 
Systems  Incorporated.  The  sensors  were  mounted  on 
the  horizontal  booms  adjacent  to  the  three-dimensional 
pitot  probes  for  correlation.  A  typical  installation 
is  shown  in  Figure  13. 

The  sensors  are  designed  for  use  with  a  constant- 
temperature  anemometer  circuit,  which  measures  the 
instantaneous  heat  transfer  rate  between  the  heated 
sensor  and  the  surrounding  air.  The  circuit  includes 
an  electrical  current  through  the  sensor  to  corre¬ 
spond  to  the  cooling  effect  of  the  flowing  air  while 
maintaining  the  sensor  at  a  constant  resistance  and 
temperature. 

These  six  anemometer  circuits  were  packaged  in  module 
form  and  were  supplied  with  DC  power  from  a  Thermo- 
Systems  module  power  supply.  In  order  to  reduce 
system  weight,  a  time-sharing  system  of  anemometer 
modules  was  used.  As  will  be  discussed  subsequently, 
three  probes  (six  channels  of  data)  on  one  side  of 
the  helicopter  (left  boom)  were  used  through  the 
anemometer  circuits  fed  through  voltage-controlled 
oscillators,  multiplexed,  and  recorded  for  30  seconds. 
Then  by  electronic  switching,  the  three  probes  on  the 
other  side  of  the  helicopter  (right  boom)  were  fed 
through  the  same  anemometer  circuits  and  recorded 
for  the  next  30  seconds. 

b.  Steady-State  Velocities 

Three-dimensional,  directional  pitot  probes,  United 
Sensor  and  Control  Corporation  Model  DA-125  modified, 
were  used  to  measure  time-averaged  velocities.  The 
probes  with  a  stem  length  of  6  inches  and  1/8-inch 
diameter  have  a  five-hole,  prism-shaped  measuring 
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head.  The  installation  of  a  three-dimensional 

probe  is  shown  in  Figure  13. 

c.  Flight  Parameters 

(i)  Altitude  -  An  absolute  pressure  transducer, 
Statham  Model  P968-154-350 ,  was  used  for 
altitude  sensing.  Installation  of  this  unit 
is  shown  in  Figure  14. 

(ii)  Airspeed  -  A  differential  pressure  transducer, 
Statham  Model  IW96TC+. 5-350,  was  used  for 
sensing  airspeed.  This  unit  was  coupled  into 
the  existing  pressure  lines  of  the  test 
vehicle  airspeed  system.  Installation  is 
shown  in  Figure  14. 

(iii)  Normal  Acceleration  -  A  Statham  temperature- 
controlled  linear  accelerometer.  Model  AJ17- 
1-350,  was  used  to  measure  the  normal  accelera¬ 
tion  of  the  test  vehicle.  The  range  of  the 
instrument  is  +1.0  g,  and  it  has  a  natural 
frequency  of  35  c.p.s.  Structural  accommoda¬ 
tion  necessitated  that  the  accelerometer  be 
mounted  on  the  aft  gear  cross  strut  (station 
116).  Location  of  the  accelerometer  is  shown 
in  Figure  11. 

(Iv)  Engine  r.p.m.  -  An  Electro  Products 

Laboratories  electromagnetic  pickup,  Model 
3015A,  was  used  to  monitor  the  engine  speed 
at  the  output  shaft.  Magnetic  field  inter¬ 
ruption  at  1/rev.  was  provided  by  an  existing 
tab  washer  on  the  shaft.  The  installation  is 
shown  in  Figure  15. 

(v)  Rotor  r.p.m.  -  The  sensor  used  for  the  measure¬ 
ment  was  the  same  type  as  that  for  measuring 
engine  r.p.m.  Field  interruption  at  3/rev.  was 
provided  by  an  arm  on  the  rotating  swash  plate. 
The  installation  is  shown  in  Figure  16. 

(vi)  Manifold  Pressure  -  An  absolute  pressure 
transducer,  Statham  Model  PA208TC-15-350, 
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was  used  to  measure  the  manifold  pressure. 
This  unit  was  coupled  into  the  existing 
manifold  pressure  line  in  the  helicopter. 

The  installation  is  shown  in  Figure  14. 

(vii)  Fuselage  Pitch  Attitude  -  An  attitude  gyro, 
Giannini  Model  3111C-1,  was  used  to  measure 
the  fuselage  pitch  attitude.  This  unit  was 
installed  on  the  pilot  side  of  the  cabin 
floor,  as  shown  in  Figure  17. 

(viii)  Longitudinal  Cyclic  Displacement  -  A  position 
indicator,  Lockheed  Model  WR8,  was  used  to 
measure  the  longitudinal  cyclic  control  input 
in  inches  of  stick  travel.  The  transducer 
was  mounted  on  the  cabin  floor  and  coupled  to 
the  long  cyclic  control  shaft. 

(ix)  Collective  Stick  Displacement  -  The  trans¬ 
ducer  used  for  this  measurement  is  the  same 
as  that  used  for  the  longitudinal  cyclic. 

The  transducer  was  mounted  on  the  pylon  shelf 
and  tied  to  the  extension  shaft  of  the 
collective  stick  as  shown  in  Figure  14. 

2.  Data  Recording  Equipment 
a.  Hot  Film  Data 


A  two-track  tape  recorder,  Viking  Model  230,  was 
utilized  to  record  the  multiplexed,  time-varying 
velocity  data  on  magnetic  tape.  Figure  18  shows  the 
recorder  installation.  As  was  previously  mentioned, 
the  six,  two-component  velocity  sensors  comprising 
twelve  data  channels  were  time-shared.  This  reduced 
the  amount  of  circuitry  required,  thereby  reducing 
the  required  size  of  the  recording  instrumentation 
package.  The  switching  circuitry,  anemometer 
modules,  and  tape  recorder  package  are  shown  in 
Figure  19. 
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b. 


Time-Averaged  Pressure  Data 


A  32-tube  inclined  manometer  bank  photo  panel  was 
used  to  depict  the  time-averaged  pressure  data  sensed 
by  the  three-dimensional  pressure  probes  positioned 
in  the  rotor  downwash.  These  data  were  recorded  by 
the  use  of  a  Graphlex  camera  equipped  to  handle 
Polaroid  film.  Figure  18  shows  the  manometer  box 
installation  in  the  helicopter. 

c.  Oscillograph 

A  recording  oscillograph,  Century  Electronics 
Corporation  Model  409E,  coupled  to  a  Century  Model 
1809  bridge  control  unit,  was  used  for  recording  the 
flight  parameters  given  on  page  48.  The  paper  speed 
used  was  from  10  in. /sec.  to  1/2  in. /sec.,  depending 
upon  the  tests.  Figure  17  shows  the  oscillograph 
installation. 

3.  Calibration  of  Instrumentation 


a.  Transducer  Bench  Calibration 


A  summary  of  the  transducer  calibration  is  presented 
in  Table  II.  The  transducer  bench  calibration  was 
performed  as  follows: 

(i)  Altitude 

Static  calibrations  were  performed  by 
utilizing  a  vacuum  pump  and  a  barometer. 

The  calibration  curve  is  shown  in  Figure  20. 

(ii)  Airspeed 

Airspeed  calibrations  were  performed  on  the 
pressure  transducer  utilizing  a  variable 
pressure  air  supply  and  a  manometer.  The 
calibration  curve  is  shown  in  Figure  21. 
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(iii)  Normal  C.G.  Accelerometer 

Static  calibration  was  obtained  by  measuring 
the  accelerometer  output  for  the  normal  and 
inverted  positions.  Zero  output  was  obtained 
by  placing  the  accelerometer  on  its  side 
normal  to  the  direction  of  excitation.  The 
calibration  is  shown  in  Figure  22. 

(iv)  Rotor  r.p.m. ,  Engine  r.p.m. 

No  calibrations  were  conducted  since  outputs 

were  read  off  the  oscillogram  as  inter¬ 
ruptions  per  unit  time. 

(v)  Manifold  Pressure 

The  manifold  pressure  transducer  was 
calibrated  using  a  vacuum  pump  and  a 
manometer.  The  calibration  is  shown  in 
Figure  23. 

(vi)  Pitch  Attitude  Gyro 

Static  calibrations  were  performed  on  an 
-  attitude  table  for  several  angles  with  a 
range  of  +16°.  The  calibration  curve  is 
shown  in  Figure  24. 

(vii)  Collective  Pitch  and  Cyclic  Stick  Displace¬ 
ment 

The  electrical  output  of  the  potentiometer 
bridge  circuit  versus  stick  motion  was 
measured  on  the  helicopter  for  Doth  control 
sticks  by  displacement  of  the  control  sticks 
from  one  extreme  to  the  other.  The  calibra¬ 
tion  curves  of  the  collective  and  cyclic 
stick  are  shown  in  Figures  25  and  26, 
respectively. 

(viii)  Hot  Film  Anemometer 


The  anemometer  system  was  bench-calibrated 
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using  the  probe  calibrator,  Thermo-Systems 
Model  1125,  for  a  velocity  ranging  from  0 
to  5000  f.p.m.  The  probe  calibrator  consists 
of  a  pressure  chamber  of  low  turbulence 
supplied  with  high-pressure  air  and  a  test 
nozzle.  The  probe  to  be  calibrated  is 
placed  at  the  nozzle,  and  velocities  are 
applied  by  adjusting  a  pressure  valve  on 
the  high-pressure  air  input.  The  velocity 
at  the  probe  (or  nozzle)  is  measured  on  an 
inclined  manometer.  The  voltage-controlled 
oscillators  (VCO's)  were  set  to  measure  center 
frequency  at  3000  f.p.m.,  in  order  to  obtain 
the  best  operation  in  the  useable  portion  of 
the  VCO's.  The  probe  output  voltage  per 
channel  was  measured  with  a  Fluke  meter  and 
plotted  versus  f.p.m.  Hot  film  calibrations 
are  given  in  Figures  27  through  32. 

(ix)  Three-Dimensional  Pitot  Sensors 


The  United  Sensor  three-dimensional  pressure 
probes  were  calibrated  by  the  manufacturer. 

A  typical  set  of  calibration  curves  is  shown 
in  Figures  33  through  35. 

b.  System  Calibration  and  Operation 

Ground  calibration  of  the  helicopter  instrumentation 
was  conducted  prior  to  each  sequence  of  flight  tests. 
All  systems  were  turned  on,  set  on  standby,  and 
allowed  a  30-minute  warm-up  period  to  stabilize  the 
transducer  circuits.  All  transducer  outputs  were 
recorded  by  the  Century  oscillograph  with  the 
exception  of  the  probes. 

Those  transducers  monitored  with  the  oscillograph 
were  resistance-calibrated  using  the  manual  triggering 
circuit  of  the  Model  1809  Bridge  control  unit.  There 
was  no  resistance  calibration  for  the  r.p.m.  circuits 
since  they  were  self -energizing.  The  anemometer 
circuits  were  calibrated  for  zero  velocity  over  the 
probes  using  manual  triggering  of  the  individual 
circuits.  The  calibrations  were  recorded  on  magnetic 
tape. 
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The  instrumentation  system  was  manually  activated 
prior  to  the  test  flights.  Once  the  helicopter  had 
become  stable  in  the  airborne  test  condition,  a 
remote  control  switch  on  the  cyclic  stick  was  pilot- 
actuated  to  operate  the  recorders.  Turning  the 
recorders  off  at  termination  of  the  test  actuated 
the  automatic  Real  sequencing  switch,  thereby  giving 
in-flight  resistance  calibration. 

D.  TEST  PROCEDURES 


The  complete  test  program  is  presented  in  Table  III.  The 
test  procedure  for  each  flight  test  condition  was  as  follows 

1.  Evaluation  of  Helicopter  Gross  Weight 

The  empty  weight  and  balance  of  the  instrumented  helicopter 
were  determined  from  actual  weight  reactions  measured  at  the 
prescribed  jack  points.  The  maximum  gross  weight  attainable 
within  the  center  of  gravity  limits  was  then  calculated.  A 
typical  weight  and  balance  statement  is  shown  in  Table  III 
for  the  balance  diagram  in  Figure  36. 

Prior  to  each  day  of  flight  tests,  and  as  determined  by  the 
cumulative  number  of  test  hours  completed,  the  fuel  load  was 
topped  off,  giving  a  known  test  gross  weight.  Having  the 
known  gross  weight,  the  total  elapsed  time,  and  the  required 
gallons  of  fuel  to  fill  the  tank  for  the  given  number  of 
hours  flown,  an  average  fuel  consumption  rate  was  calculated 
for  determination  of  test  gross  weight  per  flight.  The 
flight  test  configurations  tabulated  in  Table  IV  were 
conducted  at  a  nominal  gross  weight  of  1600  pounds.  The 
pressure  probes  and  hot  film  sensors  were  aligned  as  shown 
in  Figure  13. 

2.  Checkout  of  Instrumentation 


The  instrumentation  recorders  and  the  transducer  bridge 
circuits  were  subjected  to  a  1/2-hour  warm-up  period  to 
stabilize  the  circuitry.  Instrument  zeros  and  resistance 
calibrations  were  then  measured  and  manually  recorded.  The 
zeros  and  resistance  calibrations  were  conducted  prior  to 
each  test  flight. 
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The  instrumentation  was  activated  by  ground  personnel  prior 
to  each  flight,  and  the  instrument  recorders  were  activated 
during  flight  by  a  toggle  switch  which  was  attached  to  the 
cyclic  control  stick.  The  recording  sequence  conducted  by 
the  pilot  was  as  follows: 

a.  Activate  recorders  by  placing  switch  in  up  position. 

b.  After  a  period  of  approximately  30  seconds,  as 
advised  by  ground  control,  move  switch  to  full-down 
position  to  activate  photo  manometer  camera. 

c.  Immediately  return  switch  to  full-up  position; 
record  for  30  seconds. 

d.  Return  switch  to  off  position.  Returning  the  switch 
to  off  position  activated  a  time-delay  relay, 
allowing  the  resistance  calibration  auto-sequencing 
switch  to  complete  its  full  cycle. 

At  the  termination  of  each  test,  the  transducer  and  Real 
zeros  were  recorded  manually.  The  magnetic  tape  footage 
used  during  each  test  was  also  recorded. 

3.  Flight  Test  Procedure 

a.  Hovering  Out  of  Ground  Effect 

The  measurements  of  instantaneous  and  average  induced 
velocities  were  made  at  the  specified  flight  condi¬ 
tions,  at  three  vertical  locations  and  six  azimuthal 
positions.  The  measurements  wTere  investigated 
simultaneously  with  time-sharing,  as  discussed  in 
Section  III.  C.  Longitudinal  and  lateral  balance 
was  provided  by  testing  azimuthal  position  1A  with 
1R,  2A  with  2R,  and  3A  with  3R,  as  shown  in  Figure 
37.  The  test  rig  was  mounted  along  the  centerline 
of  the  helicopter  lateral  axis,  parallel  to  the  rotor 
shaft.  The  horizontal  booms  were  fixed  in  vertical 
position  I  and  azimuthal  positions  1A  and  1R.  The 
helicopter  then  ascended  to  the  test  altitude  of 
5000  feet.  When  the  helicopter  was  trimmed  out  at 
a  steady-state  hovering  condition,  the  data  were 
taken.  The  helicopter  then  landed  for  the  next  test. 
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While  on  the  ground,  the  horizontal  booms  were 
rotated  and  fixed  in  azimuthal  positions  2A  and  2R 
at  the  same  vertical  position  I ,  and  the  complete 
procedure  was  repeated.  These  series  were  repeated 
for  different  vertical  heights  and  azimuthal 
positions,  as  given  in  Table  IV. 

b.  Hovering  In  Ground  Effect 

An  indexed  drag  chain  attached  to  the  skid  was  used 
as  a  visual  reference  to  determine  altitude  during 
hovering  tests  while  in  ground  effect.  Thus,  the 
ground  observer,  in  audio  control  with  the  pilot,  was 
able  to  direct  his  hovering  altitude.  Hovering 
altitudes  of  zero,  3,  and  7  feet  were  used  for 
testing,  as  shown  in  Table  IV.  The  procedure  and 
data  obtained  are  similar  to  those  of  the  OGE  tests. 

c.  Forward  Speed 


The  test  procedure  for  forward  speed  was  essentially 
the  same  as  that  for  hovering.  Since  in  forward 
flight  the  rotor  wake  was  skewed  at  an  angle  to  the 
vertical  shaft  line,  the  test  rig  was  moved  to 
longitudinal  station  II  in  order  to  obtain  measure¬ 
ment  in  the  rotor  wake. 

4.  Flight.  Test  Log 

A  flight  test  data  sheet,  as  shown  in  Figure  38,  was  logged 
for  each  test  giving  the  flight  number,  weather,  configura¬ 
tion,  helicopter  parameters,  estimated  gross  weight,  and 
elapsed  time  of  the  test.  The  helicopter  parameters  were 
obtained  from  the  pilot  via  radio  communication  during 
flight.  However,  this  practice  was  discontinued  partway 
through  the  hovering  test  sequence  at  the  request  of  the 
field  control  tower.  Subsequent  data  was  received  from  the 
pilot  at  the  termination  of  each  test. 

For  each  of  the  flight  conditions  described  in  item  (3) 
above,  the  following  helicopter  flight  parameters  were 
recorded : 
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a.  altitude 

b.  airspeed 

c.  normal  acceleration 

d.  engine  r.p.m. 

e.  main  rotor  r.p.m. 

f.  manifold  pressure 

g.  fuselage  pitch  altitude 

h.  longitudinal  cyclic  stick  position 

i.  collective  pitch  stick  position 

E.  DATA  REDUCTION 
1.  Time-Varying  Velocities 

The  hot  film  anemometer  signal  data  were  obtained  from  the 
magnetic  tape  by  playing  the  multiplexed  signals  through  a 
selector  bandswitching  discriminator,  Sonex  Model  34,  and 
by  recording  the  selected  channel  output  with  a  Brush 
recorder,  Model  Mark  II. 

The  discriminator  center  frequency  settings  used  were  those 
set  during  the  bench  test  calibrations.  A  fluke  meter  was 
used  to  set  the  zero-velocity  output  levels. 

Two  Brush  recorder  speeds  were  used  in  the  data  recovery. 

All  test  data  obtained  were  reduced  at  a  paper  speed  of  2.0 
mm. /sec.;  selected  tests  were  recorded  at  50.0  mm. /sec.  to 
facilitate  frequency  and  peak  measurements  of  the  flow 
velocity.  Typical  Brush  recordings  are  shown  in  Figure  39. 
The  characteristics  of  the  flow  field  can  be  obtained  as 
follows,  referring  to  Figure  39: 

The  readings  obtained  from  channels  1  and  2  at  a  time  ti 
from  the  correlation  pulse  are  3.7  and  3.9  volts, 
respectively.  As  can  be  seen  from  the  calibration 
curves  for  this  probe,  Figure  27,  zero  velocity  is 
represented  by  a  5-volt  level;  therefore,  the  velocity 
voltages  for  channels  1  and  2  would  be  5  volts-2.5  volts 
and  5  volts-3.9  volts  or  2.5  and  1.1  volts,  respectively. 
These  voltages  represent  velocities  of  600  ft. /min.  and 
1300  ft. /min.  for  channels  1  and  2,  as  obtained  from 
Figure  27.  The  instantaneous  velocity  at  ti  is  given  by 


55 


v=yVl2+ v22 

=  n/ 6002+  I3002 
=  1430  ft/min 

2.  Steady-State  Velocities 

The  three-dimensional  probe  data  were  recovered  from  measure¬ 
ments  of  the  manometer  tube  fluid  levels  recorded  on  Polaroid 
film.  A  typical  data  photograph  is  shown  in  Figure  40.  The 
data  were  processed  by  the  following  equation  to  obtain  flow 
directions  and  average  pressure: 

Measured  values  uncorrected  for  manometer  board 
incliniation: 


P,  1.0 
P2  -0.27 
P3  0.378 
P4  0.97 
P5  0.594 
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The  yaw  angle  4  is  determined  from  the  ratio 


P2-P3 


P.-P 


0.685 


where 


p.hlh 

2 

=  0.054 


(115) 


Referring  this  ratio  to  the  typical  calibration  curve  in 
Figure  33,  the  yaw  angle  4  is  found  to  be  -12.0°.  The 
pitch  angle  9  is  determined  from  the  ratio 

P  -  P 

-—^  =  0.397  (116) 


Referring  this  ratio  to  the  calibration  curve  in  Figure  34 
and  using  the  above  angle  4  ,  the  pitch  angle  9  is  found 
to  be  -12.8°. 

The  dynamic  pressure  coefficient  q  was  determined  from  the 
ratio 


Pq-P 
P,  -P 


(117) 


where 


q  =  po  - p 


(118) 
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Having  the  values  of  6  and  4>  ,  the  ratio  for  equation  (117) 
was  found  to  be  1.23  from  Figure  35. 

The  test  dynamic  pressure  q  was  calculated  from 


q  =  AqSj  puj 


(119) 


where 


A  =  correction  factor  for  manometer  inclination  of 
30°  =  1/2 

q  =  1.162 

Sf  =  specific  gravity  of  fluid  =  .826 
=  density  of  water 

Thus, 

q  =  3.0  lb./ft.2 


The  velocity  is  given  by 


where 


P=P0 


(120) 


(121) 
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The  velocity  is  calculated  to  be 


V  =  / - — —  (3.0) 

V. 00238 

=  50  2  ft/sec  for  T=75°F 


3.  Helicopter  Parameters 

The  transducer  data  measurements  of  the  helicopter 
parameters  were  recorded  with  the  Century  oscillograph. 

The  resulting  oscillograms  were  reduced  in  the  conventional 
manner  using  the  calibration  data  in  Table  II  and  in  Figures 
20  through  26,  as  applicable,  and  corrected  for  Real 
deflections  in  a  given  test.  A  typical  oscillogram  is  shown 
in  Figure  41. 
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IV.  TEST  RESULTS  AND  CORRELATION  WITH  THEORY 


A.  HOVERING  OUT  OF  GROUND  EFFECT 
1.  Experimental  Data 

The  results  of  the  hovering  tests  out  of  ground  effect  are 
presented  in  Figures  42,  43,  and  44.  Pressure  data  as  well 
as  time-averaged  hot  film  data  are  also  included  in  these 
figures.  The  induced  velocity  components  presented  are 
nondimens ionalized_ by  the  average  normal  induced  velocity  at 
the  rotor  center,  u0  }  as  given  by  equati  n  (79). 

The  normal  induced  velocity  ratio,  Un^u°,  is  seen  to  increase 
with  radial  distance,  r/R.  The  magnitude  at  any  radial 
distance  is  seen  to  increase  with  ax! al  distance  aft  of  the 
rotor,  x/R.  Maximum  values  of  2.0  uo  were  measured.  The 
different  values  of  the  data  for  different  azimuth  positions 
are  believed  to  be  due  to  the  ambient  wind  conditions 
encountered  during  the  tests.  The  time-averaged  hot  film 
data  indicate  similar  variations  with  r/R,  but  are  of  lower 
magnitude  than  the  data  obtained  from  the  pressure  measure¬ 
ments.  In  part,  this  may  be  due  to  the  fact  that  the  hot 
film  probes  are  two-dimensional  as  compared  to  the  three- 
dimensional  pressure  probes.  The  magnitude  of  the  tangential 
induced  velocities  is  on  the  order  of  0.1.  The  critical 
effect  of  ambient  wind  conditions  limited  the  quantity  of 
useful  data.  The  negative  sign  of  the  tangential  induced 
velocities  indicates  that  the  rotor  wake  turns  in  the 
direction  of  the  propeller  rotation. 

The  radial  induced  velocity  data  show  an  increasing  trend  with 
radial  distance,  and  a  decreasing  trend  with  axial  distance^, 
from  the  rotor  plane.  Maximum  values  of  approximately  0.5  u0 
were  obtained. 

Hot  film  data  for  the  induced  velocities  are  presented  in 
Figures  45  through  57.  Each  figure  presents  the  two-channel 
data  for  three  radial  positions.  The  actual  velocities  can 
be  obtained  by  use  of  the  calibration  data  presented  in 
Figures  27  through  32  and  equations  in  Section  III.  Average 
values  of  the  normal  induced  velocities  are  indicated  in  each 
figure.  These  average  values  together  with  the  pressure- 
derived  data  are  also  plotted  on  Figures  42,  43,  and  44.  The 
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double  amplitudes  of  the  fluctuations  from  the  average  values 
are  approximately  3  ft. /sec.  No  harmonic  analyses  have  been 
made  of  the  data.  Predominant  frequencies,  however,  appear 
to  be  at  approximately  8  and  21  cycles  per  second,  which 
correspond  to  1  and  3  cycles  per  rotor  revolution. 

2.  Correlation  of  Theory  With  Test  Data 

Figure  43  presents  a  comparison  of  the  theory  discussed  in 
Section  II  with  the  test  data.  Examining  _the  calculated 
values  for  the  normal  induced  velocities,  un  ,  it  is  noted 
that  the  triangular  blade-loading  theory,  as  expected,  gives 
results  which  more  closely  agree  with  the  trend  of  the 
experimental  data.  The_ constant  blade-loading  theory  results 
in  a  constant  value  of  un  .  The  correction  for  wake 
contraction  is  obtained  by  use  of  the  data  of  Reference  9. 

For  an  axial  distance  x/R  =  0.26,  the  wake  radii  j  is  given 
as  approximately  85 °L  of  the  blade  radius.  Performing  the 
radial  correction  only,  i.e., 


no  change  in  the  velocity  results,  but  of  course,  the 
constant  value  goes  only  up  to  0.85  R.  The  use  of  the  axial 
as  well  as  the  radial  correction,  i.  e., 


UM 

N CORRECTED 
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results  in  a  slight  increase  of  the  velocity  ratio.  More 
significant  changes  are  seen  to  occur  with  the  use  of  the 
triangular  blade  loading  theory.  Without  correction  for  wake 
contraction,  the  ratio  uN//°o  is  seen  to  follow  the  test  data 
trend  of  an  increase  with  r/R  but  with  appreciably  lower 
values.  Using  the  radial  wake  correction  only  results  in 
substantially  better  correlation  with  the  test  data.  The 
values  obtained  by  using  both  x/Re  and  r/R e  in  the  calcula¬ 
tions  provide  the  best  correlation  with  the  average  of  the 
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test  data.  It  is  of  interest  to  note  here  that  although  in 
accordance  with  the  data  of  Reference  9  a  wake  radius  of  0.85 
has  been  utilized  in  the  calculations,  the  present  pressure 
data  show  that  the  wake  radius  for  x/R  =  0.26  is  equal  to  or 
larger  than  0.9  R. 

A  comparison  of  the  theory  with  test  data  shows  that  the 
calculations  performed  with  x/Re  and  r/Re  result  in  the  lowest 
values  of  the  tangential  velocity  ratio  ua^o*  The  corre¬ 
lation  with  the  limited  test  data  is  seen  to  be  fair.  It 
should  be  noted  in  this  context  that  the  theory  does  not 
account  for  the  swirl  effects  resulting  from  the  profile 
drag.  It  is  therefore  expected  that  the  test  data  give 
higher  negative  values  of  ua  than  those  computed. 

The  variation  of  the  radial  induced  velocity  ratio  with  r/R 
for  the  different  analytical  approaches  is  also  shown  in 
Figure  44.  In  comparing  the  various  analytical  approaches, 
in  particular  for  the  normal  induced  velocity  ratio,  it  is 
concluded  that  the  method  using  x/Re  and  r/Re  wake  correction 
provides  the  best  correlation  with  the  test  data  for  the 
hovering  condition.  This  method  is  also  used  in  the 
comparisons  with  the  test  data  shown  in  Figures  42  and  44. 

It  may  be  seen  that,  considering  the  scatter  of  the  experi¬ 
mental  data,  the  calculated  values  provide  a  fair  to  good 
prediction  of  the  velocities  induced  in  the  rotor  wake  of  a 
helicopter  hovering  out  of  ground  effect. 

Also  shown  in  Figures  42  and  44  are_the  digital  computer- 
derived  results  of  Reference  1  for  uN/ii0  .  It  may  be  noted 
that  the  results  of  the  much  simpler  theory  presented  herein 
compare  well  with  those  of  Reference  1. 

The  blade  position  indicator  used  during  the  testing  unfortu¬ 
nately  gave  erratic  readings,  and  hence  it  is  not  possible  to 
compare  the  calculated  and  the  measured  values  for  the  same 
blade  position.  The  average  values  and  the  fluctuation  from 
the  average  can,  however,  be  used  to  obtain  an  indication  of 
the  accuracy  of  the  theory. 

The  calculated  instantaneous  values  for  three  blade  positions 
corresponding  to  Flight  No.  5  (x/R  =0.26,  ^  =  150°)  are 

presented  in  Table  VI.  The  corresponding  hot  film  data  are 
shown  in  Figure  49a.  Comparing  the  calculated  values  of  uN/u0 
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with  the  test  data,  it  may  be  noted  that  fair  to  good  agree¬ 
ment  is  obtained  for  both  average  values  and  double  amplitude 
of  the  velocity  fluctuations  from  the  average  values. 

B.  HOVERING  IN  GROUND  EFFECT 


1.  Experimental  Data 

The  pressure  data  of  the  normal  induced  velocity  of  the  rotor 
hovering  at  three  altitudes  in  ground  effect  are  presented  in 
Figure  58.  Also  presented  are  the  corresponding  values 
obtained  for  the  out-of-ground-ef feet  condition.  The  data 
are  given  for  each  of  the  radial  and  azimuth  positions  tested. 
The  test  data  show  that,  as  expected,  the  normal  induced 
velocity  decreases  with  a  decrease  of  hovering  altitude. 

The  variation  with  radial  distance  is  not  as  significant  in 
ground  effect  as  it  is  out  of  ground  effect.  Again,  changes 
due  to  azimuth  position  are  noted.  These  changes  are 
primarily  due  to  the  prevailing  surface  winds  which  were  on 
the  order  of  0  to  12  ft. /sec.  The  radial  and  tangential 
velocities  in  the  rotor  wake  are  much  more  affected  by  the 
ambient  wind  conditions  because  of  their  relatively  smaller 
magnitudes.  The  inaccuracies  experienced  in  the  measurement 
of  these  velocity  components  therefore  did  not  warrant  the 
presentation  of  these  data. 

Hot  film  data  of  the  induced  velocities  in  ground  effect  are 
presented  in  Figures  59  to  61.  Again,  the  velocity  fluctua¬ 
tions  are  seen  to  be  relatively  small. 

2.  Correlation  of  Theory  With  Test  Data 

In  Reference  10,  an  empirical  relationship  of  the  ratio 
between  the  average  normal  induced  velocities  at  the  rotor 
plane  in  ground  effect  and  out  of  ground  effect  was  presented. 
This  relationship,  which  was  based  on  the  data  of  References 
7  and  8,  was  represented  by  a  factor  kg  =  1-0. 09e~2iH/R)  ^  ^ 

plot  of  kg  is  shown  in  Figure  62.  Also  shown  in  this  figure 
are  values  of  kg  obtained  from  the  pressure  data  of  the 
present  program. 

It  may  be  noted  that  good  agreement  is  obtained  for  the  data 
closest  to  the  rotor  plane,  x/R  =  0.22.  For  a  given  hovering 
height,  H/R,  the  normal  induced  velocity  is  seen  to  decrease 
with  axial  distance  behind  the  rotor.  This  is  consistent 
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with  the  fact  that  the  actual  wake  expands  as  the  flow 
approaches  the  ground.  For  the  same  mass  flow,  this 
expansion,  of  course,  results  in  a  reduced  velocity. 

C.  FORWARD  FLIGHT 

1 .  Experimental  Data 

The  time-averaged  values  of  normal  un  ,  longitudinal  ul  ,  and 
lateral  uy  velocities  obtained  from  the  flight  test  are 
presented  in  Figures  63  through  68.  The  data  are  presented 
for  forward  flight  speeds  of  20  m.p.h.  and  46  m.p.h.  The 
velocities  are  nondimensionalized  by  u0  ,  as  obtained  from 
equation  (79),  and  are  preserted  versus  radial  position 
measured  from  the  hub  centerline. 

These  data  were  obtained  with  the  pressure  probes  aligned 
with  the  flight  direction  for  all  azimuth  positions.  Hence, 
the  velocity  components  measured  are  parallel  to  the  rotor 
shaft,  uN }  and  in  the  longitudinal  and  lateral  directions, 
ul  and  Uy  ,  respectively.  These  velocity  measurements  are 
for  the  total  velocity  components  and  hence  include  the 
induced  velocities  as  well  as  the  components  of  the  forward 
flight  speed.  The  helicopter  flight  and  control  parameters 
associated  with  these  data  are  presented  in  Table  VII.  _The 
data  of  the  velocity  component  normal  to  the  fuselage,  uN  t 
for  ^  =  30°  and  x/R  =  0.22,  show  maximum  values  of  2.5  Uo 
and  4.0  Uo  for  forward  speeds  of  20  and  46  m.p.h.,  respec¬ 
tively.  The  corresponding  velocities  at  the  other  azimuth 
positions  are  substantia1 ly  lower.  The  fuselage  attitude 
during  these  tests  was  on  the  order  of  one  degree,  and  hence 
the  forward  speed  contribution  to  these  measurements  is 
relatively  small.  At  larger  distances  below  the  rotor,  i.e., 
for  increasing  values  of  x/R,  the  variation  of  un^uo  with 
azimuth  is  seen  to  decrease.  The  longitudinal  velocity 
measurements  are  generally  seen  to  be  on  the  order  of  the 
forward  speed,  V.  Some  reductions  from  the  value  of  V  are 
noted.  These  are  believed  to  be  primarily  due  to  the  effect 
of  flow  retardation  caused  by  the  rotor  drag. 

The  lateral  velocities,  uy  }  include  only  a  minor  component 
of  the  forward  speed.  In  fact,  if  it  is  assumed  that  rotor 
lateral  flapping  can  be  neglected,  uy  =  /u2  +  u|  « ur  .* 
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The  measured  value  of  the  ratio  Uy  ^u°  is  on  the  order  of 
0.4,  which  is  similar  to  that  measured  for  the  hovering 
conditions . 

Hot  film  data  for  the  two  forward  speeds  are  presented  in 
Figures  63  through  68.  The  average  values  of  these  data  are 
generally  in  agreement  with  the  pressure  data.  The  fluctua¬ 
tions  from  the  average  values  are  seen  to  be  smaller  than 
those  measured  for  the  hovering  flight  conditions. 

2.  Correlation  of  Theory  With  Test  Data 

Some  theoretical  values  of  un,/uo  were  computed,  and  the 
results  are  shown  in  Figure  63  for  ^  =  30  and  270°.  As 
expected,  this  theory,  which  does  not  account  for  azimuthal 
variation  of  vortex  strength,  fails  to  predict  the  rather 
large  variation  of  induced  velocities  as  measured.  The 
calculated  values  are  lower  for  ^  -  30°  and  higher  than  the 
test  data  for  if/  =  270°,  and  appear  to  give  an  average  of  the 
measured  values.  As  the  axial  distance  from  the  rotor 
increases,  the  applicability  of  the  theory  may  be  expected 
to  increase.  This  is  also  borne  out  by  the  previously 
mentioned  decrease  of  the  azimuthal  variations  for  increasing 
x/R. 

Some  computed  instantaneous  values  for  the  normal  induced 
velocities  are  presented  in  Table  VIII.  It  may  be  noted  that 
the  fluctuations  from  the  average  values  are  smaller  than 
those  computed  for  the  hovering  conditions.  This  seems  to  be 
confirmed  also  by  the  test  data  of  Figures  69  through  71. 
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uN  AXIAL  VELOCITY  COMPONENT 
ua  TANGENTIAL  VELOCITY  COMPONENT 
ur  RADIAL  VELOCITY  COMPONENT 


Figure  1.  Definitions  of  Velocity  Components  and 
Coordinates  of  Point  P  in  the  Rotor 
Wake  for  Axial  Flow  Conditions. 
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CENTRAL  VORTEX 


Figure  2.  Representation  of  the  Flow  Field  of  a 

Rotor  Blade  by  a  Vortex  System  -  Axial  Flow 


I 

I 

I 


Y  Z 


Figure  3.  Velocity  Induced  by  the  Central 
Vortex  at  Point  P. 


68 


Vcos  <fc 


Figure  6.  Rotor  in  Yawed  Flow. 
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Figure  9 
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Graphical  Representation  of  Wake  Contraction. 
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Figure  10.  Installation  of  Test  Booms  and  Support 
Structure  on  Test  Helicopter,  Side  Vie 
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Figure  12.  Data  Sensing  and  Acquisition  System. 
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Figure  13.  Velocity-Sensing  Probes  Mounted  on 
Horizontal  Boom. 
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of  Hughes-269  Helicopter. 


Figure  15.  Engine  R.P.M.  Magnetic  Pickup  Instrument  on  Output  Shaft. 
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Figure  16.  Rotor  R.P.M.  Magnetic  Pickup  Instrument  on  Swash  Plate. 


Figure  17.  Fuselage  Pitch  Gyro  Mounted  on  Cabin  Floor. 


Figure  19.  Recorder  and  Anemometer  Module  Package. 
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Figure  20.  Altitude  Calibration 


STATHAM  TRANSDUCER  PM96TC  1  .3-330 
BRID6E  VOLTAGE  =  IOV  DC 
T  =  24 .  a*c. 

HS  =  29.8  IN. 

RCAL  =  100  K 

RCAL  LOAD  =  3.22  IN.  OF  WATER 


UH9IN  33VH1 ) 

Nl  *  g  'NOI103Td3a  30V81 


87 


t 


Figure  21.  Airspeed  Calibration. 
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Figure  22.  Accelerometer  Calibration. 
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Figure  23.  Manifold  Pressure  Calibrate 
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Figure  24.  Helicopter  Fuselage  Pitch  Calibration. 
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Figure  25.  Collective  Stick  Calibration 
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Figure  26.  Longitudinal  Cyclic  Stick  Calibration. 
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Figure  28.  Hot  Film  Anemometer  Calibration,  Probe  No 
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Figure  31.  Hot  Film  Anemometer  Calibration,  Probe  No. 
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Figure  32.  Hot  Film  Anemometer  Calibration,  Probe  No. 
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Figure  33.  3-D  Pressure  Probe,  Yaw  Angle  Calibration. 
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Figure  34.  3-D  Pressure  Probe,  Pitch  Angle  Calibration. 
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Figure  35.  3-D  Pressure  Probe,  Velocity 

Head  Calibration. 


FLIGHT  TEST  REPORT 


Aircraft  Type  _  Flight  No. 


Flight  From  _  To  _  Date  _ 

Time  _  Arrive 

Total  Gross  Weight  _  _  Leave 

Flight  Conditions:  Altitude  Feet 

Cruising  Speed  MPH 

Outside  Air  Temperature  _ °F 


Rotor  r.p.m.  _ 

Engine  r.p.m.  _ 

Manifold  Pressure  _ 

Ground  Conditions:  Ambient  Pressure  _  In.  Hg. 

Ambient  Temperature  _  °F 

Wind:  Direction  _ , 

Velocity  _ 

Conditions : 

Clear  _  Rain  _ 

Dusty  _  Snow  _ 

Cloudy  _  Hail  _ 

Other  _ 

Boom  Configuration: 

Remarks : 

Figure  38.  Flight  Test  Data  Sheet 
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Figure  39.  Typical  Brush  Recordings  of  Hot  Film  Data. 


Figure  40.  Typical  Photograph  of  Steady-State  Velocity  Data. 
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Timing  Trace 
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Figure  41.  Typical  Oscillogram  of  Helicopter  Parameters. 


X/R  =0.22  06E  AZIMUTH 

R  =  0.9  R  30* 


O^oo<ao+<x* 


o  o  o  o  o 

<n  m  —  n  ro 

—  CM  c\l  rO 


108 


Figure  42.  Measured  Induced  Velocities  in  the  Rotor  Wake, 
Hovering,  x/R  =  0.22  (Time  Averaged). 
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Figure  43.  Measured  Induced  Velocities  in  the  Rotor  Wake, 
x/R  -  0.26  (Time  Averaged). 
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Figure  44.  Measured  Induced  Velocities  in  the  Rotor  Wake, 
Hovering,  x/R  =  0.49  (Time  Averaged). 
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Figure  45.  Variation  of  Instantaneous  Velocity  During  Hovering, 
OGE,  30°,  x/R  =  0.22. 
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Figure  46.  Variation  of  Instantaneous  Velocity  During  Hovering 
OGE,  t  =  150°,  x/R  =  0.22. 
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Figure  47.  Variation  o^^stantaneous  Velocity  During  Hovering, 


330 
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Figure  48.  Variation  of  Instantaneous  Velocity  During  Hovering, 
OGE,  t  =  330°,  x/R  =  0.22. 
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Figure  49a.  Variation  of  Instantaneous  Velocity  During  Hovering, 
OGE,  t  =  150°,  x/R  =  0.26. 
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Figure  51.  Variation  of  Instantaneous  Velocity  During  Hovering, 
OGE,  t  =  330°,  x/R  =  0.26. 
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Figure  52.  Variation  of  Instantaneous  Velocity  During  Hovering 
OGE,  t  =  30° ,  x/R  =  0.49. 
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Figure  53.  Variation  of  Instantaneous  Velocity  During  Hovering, 
OGE,  'f'  =  90°,  x/R  =  0.49. 
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Figure  55.  Variation  of  Instantaneous  Velocity  During  Hovering, 
OGE,  \f/  =  210°,  x/R  =  0.49. 
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Figure  57.  Variation  of  Instantaneous  Velocity  During  Hovering, 
OGE,'l'  =  330°,  x/R  =  0.49. 
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Figure  58.  Measured  Normal  Induced  Velocities  In-Ground  Effect. 
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Variation  of  Instantaneous  Velocity  During  Hovering, 
In-Ground  Effect,  h/R  =0.53. 
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Figure  62.  Effect  of  Ground  on  Rotor  Induced  Velocities. 
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Figure  63.  Variation  of  Normal  Velocity  With 

Azimuth  and  Distance  Below  the  Rotor, 
V-20  m.p.h. 
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Figure  64.  Variation  of  Longitudinal  Velocity  With 
Azimuth  and  Distance  Below  the  Rotor, 
V=20  m.p.h.  (Pressure  Data). 
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Figure  65.  Variation  of  Lateral  Velocity  With  Azimuth 
and  Distance  Below  the  Rotor,  V=20  m.p.h. 
(Pressure  Data). 
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Figure  66.  Variation  of  Normal  Velocity  With  Azimuth 
and  Distance  Below  the  Rotor,  V=46  m.p.h. 
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Figure  67.  Variation  of  Longitudinal  Velocity  With 
Azimuth  and  Distance  Below  the  Rotor, 
V=46  m.p.h.  (Pressure  Data). 
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Figure  68.  Variation  of  Lateral  Velocity  With  Azimuth 
and  Distance  Below  the  Rotor,  V-46  ra.p.h. 
(Pressure  Data). 
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Figure  69.  Variation  of  Instantaneous  Velocity  During  Forward 
Flight,  x/R  =0.26,  t  =  330°. 
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Figure  71.  Variation  of  Instantaneous  Velocity  During  Forward 
Flight,  x/R  =  0.26,  r/R  =  0.29. 


TABLE  I 

DIMENSIONS  AND  CHARACTERISTICS 

OF  TEST  HELICOPTER 

GENERAL: 

Gross  Weight 

1,600  lbs. 

Useful  Load 

590  lbs. 

Overall  Length 

28  ft.  5  in. 

Overall  Height 

8  ft.  2  in. 

Rotor  Clearance  (Ground  to  tip, 

rotor  static) 

7  ft.  2  in. 

C.G.  Limit  (Total  Travel) 

5  in. 

MAIN  ROTOR  DATA: 

Rotor  Diameter 

25.29  ft. 

Number  of  Blades 

3 

Disc  Area 

502  ft. 

Blade  Chord 

9.562  ft. 

Blade  Solidity 

0.043 

Blade  Airfoil 

NACA  0015 

Tip  Speed 

596  ft. /sec. 

Rotor  r.p.m. 

450  r.p.m. 

Disc  Loading 

3.18  lbs. /ft.2 

Rotor  System 

Single  (Articulated) 

POWER  PLANT: 

Engine 

Lycoming 

Power  (Takeoff) 

180  HP 

Power  (Normal) 

160  HP 

PERFORMANCE : 

Maximum  Forward  Speed 

75  kn. 

Endurance  (2700  Engine  r.p.m.) 

2  hrs.  @  65  kn. 

Hover  Ceiling,  STOL  Day  OGE 

4,000  ft. 

Hover  Ceiling,  110°F  Day  OGE 

1,200  ft. 

\ 
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CALIBRATION  DATA 
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TABLE  II  (Continued) 


*  Hot  film  probe  correlation  trace  on  oscillograph  for  probes  recording  on  left 
side  of  helicopter.  Trace  off  paper  for  right-siae  recording. 


TABLE  III 


TYPICAL  WEIGHT  AND  BALANCE 


Aircraft  Configuration:  22  November  1966 

Instrumentation  and  related  structure  installed 
Doors  removed 

Cabin  trim  removed  except  pilot  seat  and  back 
Zero  fuel 

Full  oil;  engine  and  transmissions 
Skid  wheels  stowed  forward 


Most  Aft  C.G. 

, ,  Booms  at 

Station  137. 

0 

SCALE 

WEIGHT  (#) 

ARM  (in.) 

MOMENT  x  10" 3 

Left  Main 

680.0 

96.0 

65.280 

Right  Main 

597.0 

96.0 

57.312 

Tail 

28.0 

267.0 

7.476 

Empty  Weight 

1305.0 

99.6 

130.068 

Pilot 

170.0 

83.9 

14.263 

Fuel  (25  Gal. 

)  150.0 

107.0 

16.050 

Gross  Weight 

1625.0 

98.7 

160.381 

B.  Most  Forward  C.G.,  Booms  at  Station  100.0 


Left  Main 

678.0 

96.0 

65.088 

Right  Main 

608.0 

96.0 

58.368 

Tail 

19.0 

267.0 

5.073 

Empty  Weight 

1305.0 

98.5 

128.529 

Pilot 

170.0 

83.9 

14.263 

Zero  Fuel 

-- 

Gross  Weight 

1475.0 

96.8 

142.792 

The  instrumented  booms  positioned  in  azimuth  are  fixed  at  180° 
relative  to  each  other  for  all  test  configurations.  Thus  the 
composite  C.G.  of  the  two  booms  is  not  affected  by  changes  in 
azimuth  positions. 

Hughes  Helicopter  Model  269A 
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TABLE  IV 


FLIGHT  TEST  SCHEDULE 


Flight 

Velocity 

Boom* 

Horizontal 

Station 

Position* 

Vertical 

Azimuth 

Position 

Purpose 

I 

OGE 

Hover 

I 

1 

1A,  1R 

Azimuthal 

Vel . 

Var. 

2 

2A,  2R 

3 

3A,  3R 

4 

2 

1A,  1R 

5 

2A,  2R 

6 

3A,  3R 

7 

3 

1A,  1R 

8 

3A,  3R 

9 

OGE 

Hover 

3A ,  3R 

10 

0-Foot 

3 

1A,  1R 

Azimuthal 

Hover 

Vel. 

Var. 

in  Ground 

Effect 

11 

3A,  3R 

12 

3-Foot 

1A,  1R 

Hover 

13 

3A,  3R 

14 

7-Foot 

1A,  1R 

Hover 

15 

I 

3 

3A,  3R 

1 

TABLE  IV  (Continued) 


Flight 

Flight 

Velocity 

Boom* 

Horizontal 

Station 

Position* 

Vertical 

Azimuth 

Position 

Purpose 

16 

20  Knots 

I 

1 

1A,  1R 

Azimuthal 

Vel . 

Var. 

in  Fwd. 

Flight 

17 

] 

[ 

2A 

II 

2R 

18 

I 

3R 

II 

3A 

19 

20  Knots 

I 

2 

1A,  1R 

20 

I 

2A 

II 

2R 

21 

I 

3R 

II 

3A 

22 

[ 

3 

1A,  1R 

23 

r 

2A 

II 

2R 

24 

[ 

3R 

II 

3A 

25 

40  Knots 

[ 

1 

1A,  1R 

26 

[ 

2A 

II 

2R 

27 

I 

3R 

II 

3A 

28 


I 


2 


1A,  1R 


TABLE  IV  (Concluded) 

Boom* 

Flight 

Horizontal 

Position* 

Azimuth 

Flight 

Velocity 

Station 

Vertical 

Position 

Purpose 

29 

I 

2A 

Vel. 

II 

2R 

Var.  in 

Fwd. 

30 

I 

3R 

Flight 

II 

3A 

31 

I 

3 

1A,  1R 

32 

40  Knots 

I 

2A 

II 

2R 

33 

I 

3R 

II 

3A 

*See  Ta 

ble  V  and 

Figure  37  for  definitions. 
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TABLE  V 
TEST  NOTATION 

Longitudinal  Position  of  Boom 
I  Station  100.0  (see  Figure  36) 
Vertical  Position  of  Boom 
I  Station  100.0 


1.  W.L.  67.0 

2.  W.L.  60.5 

3.  W.L.  25.0 


Sensor  Coordinates 

Sensor 

PR 

1 

PR 

2 

PR 

3 

Description 

iebs3 

B.L. 

QS1 

B.L. 

B.L. 

Rt.  Boom: 

— 

3A 

Aft 

172.4 

-50.2 

186.6 

-54.3 

200.8 

-58.4 

LA 

Outboard 

100.0 

-105.0 

100.0 

-120.0 

100.0 

-135.0 

2A 

Forward 

27.6 

-50.2 

13.4 

-54.3 

-0.8 

-58.4 

PL 

1 

PL 

2 

PL 

Sta. 

B.L. 

Sta. 

B.L. 

Sta. 

B.L. 

Left  Boom: 

2R 

Aft 

172.4 

50.2 

186.6 

54.3 

200.8 

58.4 

1R 

Outboard 

100.0 

105.0 

100.0 

120.0 

100.0 

135.0 

3R 

Forward 

27.6 

50.2 

13.4 

54.3 

-0.8 

58.4 
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TABLE  V 


(Continued  ) 
Longitudinal  Position  of  Boom 
II  Station  137.0 
Vertical  Position  of  Boom 
II  Station  137.0 

1.  W.L.  65.0 

2.  W.L.  62.0 

3.  W.L.  25.0 


Sensor 

Sensor 

Coordinates 

PR 

1 

PR 

2 

PR 

3 

Position 

Description 

Sta. 

emu 

Sta. 

B.L. 

B.L. 

Rt.  Boom: 

3A 

Aft 

204.0 

-62.6 

217.6 

-68.8 

230.8 

-75.0 

1A 

Outboard 

137.0 

-105.0 

137.0 

-120.0 

137.0 

-135.0 

2A 

Forward 

62.0 

-36.0 

47.2 

-37.2 

32.4 

-38.4 

PL 

1 

PL 

l 

PL 

3 

0W 

Baa 

B.L. 

H3M 

B.L. 

Left  Boom: 

2R 

Aft 

204.0 

62.6 

217.6 

68.8 

230.8 

75.0 

1R 

Outboard 

137.0 

105.0 

137.0 

120.0 

137.0 

135.0 

3R 

Forward 

62.0 

36.0 

47.2 

37.2 

32.4 

38.4 

TABLE  VI 


CALCULATED  VELOCITIES 


U°  -  25.86  ft. /sec. 

Flight  Condition  #5  (x/R  =  0.26,  ^  -  150°) 


X' 

re 

Blade 

Azimuth 

Position 

UN/U0 

ua/!io 

ur=ur 

0.306 

0.682 

A 

li 

-0.0665 

B 

1 :  1 

-0.00911 

0.19 

c 

1.31 

-0.0149 

0.788 

A 

1.67 

-0.0650 

B 

1.53 

-0.00837 

0.27 

C 

1.52 

-0.0104 

0.906 

A 

1.98 

-0.0510 

B 

1.78 

-0.00715 

0.36 

C 

1.74 

-0.00715 
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V. 


CONCLUSIONS 


The  results  of  the  work  performed  in  this  program  lead  to 

the  following  conclusions: 

1.  The  feasibility  of  obtaining  full-scale  induced  velocity 
data  by  use  of  the  measuring  techniques  utilized  was 
established. 

2.  Hovering  data  show  that  the  normal  induced  velocity 
increased  with  increasing  radial  distance.  Radial 
induced  velocities  are  on  the  order  of  one-half  the 
average  normal  induced  velocity,  but  tangential  induced 
velocities  are  less  than  ten  percent  of  this  value. 

3.  Test  data  obtained  for  hovering  while  in  ground  effect 
show  that  the  induced  velocity  reductions  are  similar  to 
those  obtained  by  earlier  researchers. 

4.  Significant  azimuth  variations  of  normal  induced 
velocities  were  measured  in  forward  flight. 

5.  Theoretically  calculated  values  of  the  induced  velocities 
show  some  correlation  with  the  hovering  test  data.  The 
theory  does  not  provide  good  correlation  with  the  forward 
speed  test  data  in  the  proximity  of  the  rotor. 

6.  Calculations  indicate  that  the  effect  of  the  fuselage  on 
the  induced  velocities  is  negligible  for  the  Model  269-A 
helicopter. 

7.  Experimental  error  for  time-averaged  velocity  values 
ranged  from  1  percent  to  5  percent. 

8.  Experimental  error  for  instantaneous  velocities  ranged 
from  3  percent  to  10  percent. 
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